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NOMENCLATURE 


Symbol 

Description 

A/D 

Analog to Digital 

AMD 

Arithmetic mean diameter 

CPM 

Continuous pulse mode 

D( 10) 

Arithmetic mean diameter 

D(20) 

Area mean diameter 

D(30) 

Volume mean diameter 

D 32) 

Sauter mean diameter 

D d 

Drop diameter 

D fc 

Drop diameter at background 

f 

Disturbance frequency 

GL 

Gray level 

A 

Wavelength 

MAGL 

Measured average gray level 

PBG 

Particle boundary gradient 

<t> 

Relative phase shift associated with P/DPA signals 

PSP 

Particle sizing program 

q 

Liquid flow-rate 

SD 

Standard deviation 

SMD 

Sauter mean diameter 

SPM 

Single pulse mode 

T 

Image threshold 

Tfc 

Image threshold just above background 

P 

Droplet velocity vector 
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Section 1 


INTRODUCTION 


Spray characterization is essential in many technologies. Improved cloud simulation for icing stud- 
ies, increased efficiency for combustion technology, and design optimization of applicator nozzles 
for industry and agriculture are only a few areas which benefit from accurate spray measurements. 
The lack of a universally accepted calibration/verification standard and operating characteristics of 
sizing instrumentation has left the questions of accuracy and repeatability in spray measurements 
unanswered. Recently, various groups (e.g., ASTM Subcommittee E29.04 on Characterization of 
Liquid Particles, 1986 Droplet Technology Workshop, etc.) have addressed the question of accu- 
racy and calibration in drop-size instrumentation, however no agreement has been reached with 
regard to methods or apparatus for standardizing drop-size measurement instruments [1]. The 
following work involves the evaluation of two instruments based on different drop-sizing techniques 
in side-by-side benchmark tests under identical operating conditions. 

The non-intrusive nature of laser/optical techniques have shown the most promise in spray char- 
acterization. Of the three major types of laser/optical techniques (i.e., imaging, doppler anemome- 
try, and laser-diffraction), the laser-diffraction method is most widely used, and probably the best 
known system is the Malvern instrument [2]. Doppler anemometry, however, is receiving more 
attention due to the recent development of Aerometric’s P/DPA, which has an increased sizing 
range (35:1) [3,4], in comparison to the (10:1) range for visibility dependent Doppler anemometers 
[5]. With the use of real-time digital image processing to perform focus discrimination without 
correction, the University of Nebraska - Lincoln (UNL) laser imaging system [6-10] has shown the 
capability for true volumetric analysis. Previously, imaging systems, e.g., Weiss et al. [11], and oth- 
ers, have used depth of field corrections based on the maximum measured drop-size to “back-out” 
the number of smaller particles in a normalized volume. Processing time can be saved using this 
method, however the assumptions may lead to errors in obtaining accurate size characteristics. The 
above techniques vary in several areas; 1) sampling method (e.g., spatial vs. temporal), 2) probe 
volume (e.g., line of sight averaging, crossed beams, vs. focus volume), 3) instrument drop-size 
range and resolution, and 4) calibration and/or verification (e.g., reticles, monodisperse droplets, 
or polydispersions). Similarities shared by the imaging technique and the laser-diffraction method 
are that both are spatial sampling methods which allows for similar calibration (i.e., calibration 
reticle [7,12]). The similarity in probe volume of Doppler anemometers and imaging systems al- 
low for verification and comparison with minimal correction. In this work, a P/DPA and a laser 
imaging system were evaluated by concurrently performing a set of baseline benchmark tests. 

According to Tishkoff [13], chairman of ASTM Subcommittee E29.04 on Characterization of 
Liquid Particles, the four major areas of concern in spray characterization are instrumentation, 
sampling, data processing, and terminology. In the following work, the emphasis of the evaluation 
was placed on instrumentation (i.e., the setup and operation of the P/DPA, a temporal sampling 
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instrument in ideal conditions, and the UNL laser imaging system, a true spatial sampling in- 
strument). The difference in data acquisition or sampling method was minimized by overlapping 
the probe volumes of the two systems [14] and analyzing a spray under steady-state conditions 
(i.e., spray characteristics remain constant with respect to time). Data processing and terminology 
of the two systems closely follow the standard practices established by ASTM [15]. Taking into 
account the above criteria, the comparison of the P/DPA and the UNL laser imaging system was 
accomplished with minimal reduction of drop-size data. 

The comparison of the P/DPA and the UNL laser imaging system is discussed in the following 
order; 1) experimental apparatus including the droplet sizing instruments, 2) procedure and op- 
erating conditions for the benchmark tests, 3) results obtained from the benchmark tests, and 4) 
conclusions as to the operation, data representation, and comparability of the two instruments. 
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Section 2 


EXPERIMENTAL APPARATUS AND PROCEDURE 


The apparatus, used in the benchmark tests, consisted of a P/DPA [3,4], a laser imaging/video 
processing system (LI/VPS) [6-10], a MOD-1 nozzle [16], air and water supply systems (AWSS), 
and the measurement instrumentation used to monitor the operating conditions of the nozzle. 
Verification tests were performed using a Berglund-Liu vibrating orifice aerosol generator (VOAG) 
[17,18]. Operating conditions of the tested apparatus and the setup parameters for the sizing 
instruments are detailed. 

2.1 P/DPA 

Phase/Doppler Particle Analyzer theory and operation are described by Bachalo et al. in several 
references [3,4], therefore, only a brief description of the P/DPA components and operation follows. 
Setup features specific to this research are detailed with special attention given to the selection of 
appropriate photo- multiplier tube (PMT) gain voltage. 

The P/DPA is a crossed beam laser Doppler anemometer (Fig. 2.1). The P/DPA transmitter 
utilizes a 10 mW He-Ne laser. The transmitter beam is split and the resulting beams are focused to 
a point by a convex lens. The Doppler fringes, formed at the crossed beam intersection, are relayed 
to the P/DPA receiver by the refracted light from a droplet passing through the crossed beam 
intersection. The P/DPA receiver uses a pair of convex lens to collect and focus the Doppler fringes 
from the passing droplet onto three PMTs, aligned parallel to the droplet’s velocity vector (u). The 
PMT voltages are filtered and amplified to remove the pedestal component of the burst and increase 
the differentiation of Doppler frequencies in the signal (Fig. 2.2). Particle size measurements are 
determined from the phase shift in the filtered Doppler signal. 

Velocity measurements are taken identically to the laser Doppler velocimeter, but the P/DPA 
is very distinct in its method of particle size measurement. Bachalo et al. [4] have shown droplet 
size (Dj) to be dependent on the relative phase shift (<f>) associated with a Doppler signal incident 
on two adjacent PMTs. 

With the operating conditions of the VOAG and the MOD-1 nozzle varying, the P/DPA also 
required adjustment in operating parameters. The following is a brief summary of the P/DPA setup 
parameters (Fig. 2.3). Parameters (A) and (B) are specified for the transmitter laser supplied by 
the manufacturer, and do not require adjustment. Hardware parameters of the P/DPA fixed for 
the duration of this work, specified according to reference [19], were; (E) the focal length of the 
transmitter lens used, was 495 //m for a measurable size range of 1 to 300 micrometers (/im), (F) 
the receiver was positioned 30° off the forward axis of the transmitter for sizing water droplets, (G) 
the refractive index was set for water, and (T) the Direct Memory Access (DMA), which allows for 
the storage of approximately 16,000 concurrent raw PMT signals for processing, was switched off 
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Figure 2.1: Phase Doppler/Particle Analyzer 
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Figure 2.3: P/DPA Setup Page 


to facilitate the comparison with the LI/VPS. For this research, the beam separation, parameter 
(D), was alternated between 25 and 12.5 mm for the different spray size distributions generated (i.e., 
the beam separation and the transmitter lens’ focal length specify the fringe spacing and number 
in the probe volume which, in turn, specifies a range of allowable drop-sizes to be measured). 
Other parameters, such as; (N) and (M) the high pass filter setting, (L) PMT voltage, (J) size, 
and (Q) velocity ranges are set according to the specific operating conditions droplet density, size 
distribution, etc.) of the VOAG or MOD-1 nozzle. The high pass filter allows only those Doppler 
signals with a frequency above a preset limit to pass on for further processing. The high pass 
filter setting is dependent on the average droplet velocity, and can be set by studying the count vs. 
velocity distribution. The selection of a high pass filter can be fine-tuned by using an oscilloscope 
to monitor the filtered PMTs for uniform signals with minimal distortion. The previous parameters 

are discussed in detail in the P/DPA operating manual. 

The PMT gain voltage was to be set at a point just prior to PMT saturation. The above was 
accomplished by studying the saturation lights connected to each PMT. The saturation lights were 
to flash intermittently 50% of the time which implied approximately 1% saturation. Following 
the above procedure in performing an analysis on a high density spray, an inordinate number of 
large drops showed up in the analysis (Fig. 2.4). The large drops were determined to be false by 
concurrent studies by the LI/VPS and previous studies by NASA on the tested nozzle. According 
to Bachalo [20], the false drops were reflections or echoes in the PMTs caused by the high density 
of the spray, therefore, the PMT voltage should be set by stepping through the PMT voltage range 
(i.e., approximately 275 to 475 volts), and studying the number vs. size distribution for a point 
where little change occurs in the distribution shape (Fig. 2.5). 

2.2 Laser Imaging/Video Processing System 

The basic architecture of the LI/VPS has been described in detail by Ahlers and Alexander [8,9]. 
Ahlers [7] performed an analysis on static particles (e.g., polystyrene microspheres) situated in the 
plane of focus of the imaging optics. Further work by Wiles [10] described a technique for focus 
classification without depth of field corrections. The implementation of a particle sizing system 
capable of performing analysis on aerosol sprays has been the focus of the current research program. 
The following discussion is divided into sections covering: 1) components and operation, 2) drop 
sizing method, 3) calibration technique to minimize uncertainty due to camera tube non-linearities, 
4) focus criteria, 5) modifications for dynamic measurements, and 6) software updates. 

2.2.1 Components 

The LI/ VPS is divided into two subsystems, a laser imaging device and a video processor . The laser 
imaging device (Fig. 2.6) components are: a COHU camera system (control unit and camera), a 
Laser Energy Inc. (LEI) laser system (power supply unit and laser), a Laser Holography Inc. (LHI) 
control system (sync circuit and laser control unit), the imaging optics, a Panasonic NV-8950 
or RCA VET650 VCR, a Panasonic TQ-2023 (A) laser/optical memory disk recorder (LDR), a 
Panasonic WJ-180 time/date generator, a Sony Trinitron monitor, a Sanyo monitor, and a back- 
up Molectron UV Series II Model UV12 (MUV12) N 2 laser. The video processor (Fig. 2.7) 
consists of a Recognition Concepts Inc. (RCI) Trapix 55/32 real-time image processor, a PDP 
11/73 computer for control, and the processing software. A LSI- 11/ 03 computer is also available 
for utility processing. 
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Figure 2.4: Reflections Caused by High Density Spray 



Figure 2.5: Drop Distribution Behavior with increasing PMT Voltage 
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Figure 2.6: LI/VPS Laser Imaging Device 
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Figure 2.7; LI/VPS Video Processor Schematic 
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The baseline sync of the laser imaging system originates with the camera control unit (CCU). The 
CCU, operating on 60 Hz (line) cycle, drives the camera at video rates (i.e., one field every 16.67 
milliseconds (ms) or one complete frame every 33.33 ms). The laser sync circuit (LSC); 1) receives 
the CCU triggering pulse, 2) uses the CCU trigger to generate a sync pulse for the laser, 3) sets 
the laser in sync with the camera process, and 4) sets the pulse rate of the laser to multiplies of 60 
Hz (e.g., 30, 15, etc.), or allows the operator to pulse the laser manually or by computer control. 

The LHI laser control unit has variable power settings with an internal sync generator. The 
LEI laser system consists of a Model N2-50 power supply and pulsed laser (A = 337 nm). The 
original system was operable within a range of 2-20 kW pulsed power and has been upgraded to 
40 kW. By changing the mirrors in the laser tube, the pulse duration of the laser can be varied 
from either 3 nanoseconds (ns) or 10 ns. A second N 2 laser (MUV12) also contains its own internal 
sync generator, but the power cannot be varied. The MUV12 (laser and vacuum pump) has a peak 
power output of 250 kW and is limited to a pulse duration of 10 ns. 

With the laser system in sync with the camera system, the object field is transferred to the 
camera by the imaging optics. A piano convex lens magnifies the object field before transferring 
the object field to the camera tube. System capabilities include a 500X and 1000X lens (i.e., 500X 
implies 800 by 800 micrometer (/zm) field of view, and 1000X implies 400 by 400 /im field of view ) 
for measurement. The video signal is than routed to a VCR where the images can be recorded 
for later viewing as a visual aid, or the images can be sent to the digital image processor. Other 
available options to the system are the use of the Panasonic time/date generator which overlays 
the time, date, and optional stopwatch capabilities on the analog video signal; and the availability 
of the Panasonic TQ-2023F LDR to store video frames which can provide for fast retrieval time 

without the tape positioning problems associated with a VCR. 

The user interfaces with the LI/VPS at the PDP 11/73 console. Through the processing 
software, the user instructs the Trapix 55/32 to perform various logical and arithmetic operations 
on the images supplied by the laser imaging system. The Trapix 55/32 image processor has one 
megabyte of image memory which gives the processor available space to store four concurrent video 
frames. The PDP 11/73 computer controls the Trapix 55/32 through a parallel interface with a 
sub-library of control subroutines. The LSI-11/03 computer is also available for utility processing. 


2.2.2 Sizing Method: Segmentation 

The original software package developed by Ahlers [7] uses a technique called segmentation. The 
segmentation technique was adopted because sequential line by line processing is inherent to the 
camera system. The camera outputs a standard RS-170 composite video signal. The video signal 
is composed of 525 scan lines with interlace (i.e., odd and even scan lines interwoven into one 
complete frame). The segmentation technique uses the pattern recognition of the system (i.e., the 
conversion of the analog video signal into discrete pixels with specific intensity level and position) 
to analyze particles. 

The premise of segmentation implies that discrete line segments, which lie adjacent to one 
another, can be summed into discrete two-dimensional objects. With the particles appearing as 
black disks on a white background in the digitized frame, the segmentation method finds the pixels 
upon whicli the particles reside and joins them into line segments (one pixel wide) in the line by 
line processing. The software matches the segments of the previous line to the current line until 
the objects are completely specified (Fig. 2.8(a)). 
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Figure 2.8: LI/VPS Particle Representation 
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The Analog- to Digital conversion is performed by the Trapix 55/32. The analog signal (i.e., video 
frame) is converted to a 512x512 array with array elements (i.e., pixels) that have eight bit pre- 
cision (i.e., 256 grey levels). Ahlers showed the optimum threshold (T) was at a gray level of 
approximately 90 [7]. Figures 2.8(b) and 2.8(c), show the digitized particle before and after the 
thresholding process has been performed, respectively. After-which, with the subroutine, FINDTR, 
developed by Ahlers [7], the processor is able to find the transition which occurs at the 90 T. With 
the two transition points of a segment found, the program processes the remainder of the line until 
all segments are found. The above procedure is the basis for segmentation with program execution 
continuing in a line by line order. 

2.2.3 Calibration 

Previous work on the LI/VPS has included sections on calibration [7,10]. The initial work by 
Ahlers determined the qualifiers for calibration and specified an initial set of magnification cor- 
rection factors (MCF). MCF qualifiers were the micron per pixel correction, the correction for 
non-linearities in the camera tube and the optimum value for the threshold of the image for sizing 
particles. The camera non-linearities initially were assumed to be dependent only on the x pixel 
location, this assumption required; 


MCF = f(x). (2.1) 

Further work by Wiles showed improved accuracy by specifying MCFs with x and y dependence; 


MCF = /(x,y). (2-2) 

In Ahlers’ work, MCFs were determined by fitting experimental data points (i.e., x position, 
MCF) to the appropriate curve (i.e., straight line, exponential, etc.), whereas with Wiles work, 
the MCFs as functions of x and y pixel position were found intuitively. In this researcher s work, 
calibration of the system became necessary after the COHU camera tube had to be replaced due 
to loss of sensitivity. Because the two-dimensional MCFs determined by Wiles were intuitive and 
specific to the replaced camera tube, a new method, which could be easily repeated, had to be 
deduced for determining the MCFs. Experimental data was discretized into 50 pixel intervals (Fig. 
2.9), whereby the MCF was implied to be constant with respect to the x position in each interval, 


MCF = < 


f i(y), 
/2(y), 
/3(y), 
/4(y), 
/ 5(y), 
/6(y)> 
/7(y)> 
/8(y)> 


50 < x < 100 
100 < X < 150 
150 < X < 200 
200 < x < 250 
250 < x < 300 
300 < x < 350 
350 < * < 400 
400 < x < 450 
for 50 < y < 450. 


(2.3) 


The above functions could than be found by curve- fitting the data (y position, MCF) specific to 
each interval. The following discussion is a description of the calibration method and procedure 
used. 

The calibration method uses a calibration reticle (i.e., opaque disks in the form of thin metal 
films deposited on glass substrate) [12]. The configuration and particle size variation of the specific 
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reticle (Model #RR-50-3.0-0.08-102-CF-114) used in calibration are shown in Fig. 2.10 and Table 
2.1. The range in diameter of the reticle particles is 5.29 /zm to 92.75 /zm. The calibration reticle is 
well suited for the LI/VPS because it can easily be positioned in the plane of focus of the imaging 
optics, eliminating the need for depth of field correction. 

The calibration procedure uses a revised version of the Particle Sizing Program (PSP) developed 
by Ahlers [7]. The modified PSP is setup to collect data (i.e., particle position, x and y pixel 
diameters, etc.) for a prescribed opaque disk from the calibration reticle. With the calibration 
reticle in the focal plane of the imaging optics, the calibration program is started. The calibration 
reticle is then positioned randomly throughout plane of focus with the program storing the data 
simultaneously. With the known diameter, the MCFs are found by Equation (2.4); 


MCF = True diameter (/xm) — (2.4) 

Measured diameter (pixels) 

The calculated MCF is then specified according to the particle’s center position. The data is 
then sorted into the perspective 50 vertical pixel intervals, and then each set of data (i.e., y pixel 
position, x MCF, and y MCF) is sorted according to y pixel position. With the correction factors 
specified as dependent variables of the y pixel position, the data can be set to the best fit curve. 
Figure 2.11 is a flow diagram of the aforementioned procedure. The above procedure was carried 
out for the 500X and the 1000X lens. The use of a different particle from the calibration reticle being 
the only change in the procedure. Because the MCFs are determined in the procedure as average 
values over the total diameter of the particle, the appropriate particle had to be chosen to avoid 
excessive overlapping of calibration intervals. Also, to avoid the edge effect (i.e., pixel elements 
being discrete implies pixels can be on or off depending on the position of the true particle s edge), 
the largest available particle should be chosen. 

Preliminary work showed that approximate MCF for the 500X lens was 2.1 /zm/pixel, and con- 
versely, 0.98 /zm/pixel for the 1000X lens. As implied above for the 50 pixel intervals, a calibration 
particle diameter of 25 pixels would minimize interval overlap and edge effects. Therefore, for 500X 
lens, the #16 particle (i.e., 52.5 /zm) was used, and conversely, for the 1000X lens, the #7 particle 
(i.e., 23.90 f/m) was used. The results of the above procedure and a comparison of previous system 
calibrations with the present calibration is presented in Section 3.1. 


2.2.4 Focus Method 

Ahlers [7] performed work using polystyrene micro-spheres restrained between two glass micro- 
scope slides positioned in the plane of focus of the imaging optics. The above tests verified the 
methodology and calibration of the LI/VPS. As with most complex systems, development occurs 
in stages, therefore Ahlers constructed a particle sizing system which performed analysis on static 
and semi-static particles in the focal plane of the imaging optics with good accuracy. Wiles [10], 
in the next stage in the development of the LI/ VPS, defined a method of focus classification (i.e., 
particles unaffected by diffraction light scatter). As Fig. 2.12 shows, with a diffraction limited 
system, particle focus is dependent on the particle’s boundary gradient and it’s relative intensity 
as compared to background. Because of the 8-bit precision of the video processor, the particle s 
intensity level with respect to background could be used as a viable criteria for focus. The parti- 
cle’s boundary gradient (PBG) was used as a secondary test because it rejects large out of focus 
particles which appear as small particles in focus by the particle intensity level test [10]. 

With the 256 grey level resolution and the processing capabilities of the video processor, the 
focus parameters are determined. The particle’s intensity level or measured average grey level 
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Figure 2.10: Calibration Reticle 
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Table 2.1: Specification Sheet for Calibration Reticle 
CALIBRATION RETICLE : RR-50-3.Q-Q,08- I02-CF - #11_4_ 
FINAL DATA SHEET 1 


DIAMETER 


AREA 

VOLUME 

(ftm) 1 2 

NUMBER 

FRACTION 

FRACTION 


1 

5.29 

2898 

0.015 

0.002 

2 

6.81 

776 

0.006 

0.001 

3 

8.98 

895 

0.013 

0.003 

4 

11.93 

1171 

0.030 

0.009 

5 

17.20 

1009 

0.054 

0.023 

6 

21.33 

642 

0.053 

0.028 

7 

23.90 

456 

0.047 

0.028 

8 

26.71 

505 

0.065 

0.043 

9 

31.11 

396 

0.069 

0.054 

10 

34.17 

280 

0.059 

0.050 

11 

37.07 

306 

0.076 

0.070 

12 

40.47 

240 

0.071 

0.072 

13 

42.71 

207 

0.068 

0.073 

14 

47.37 

160 

0.065 

0.077 

15 

50.39 

106 

0.049 

0.061 

16 

52.50 

109 

0.054 

0.071 

17 

56.23 

96 

0.055 

0.077 

18 

60.70 

88 

0.058 

0.089 

19 

67.04 

58 

0.047 

0.079 

20 

73.48 

27 

0.026 

0.048 

21 

80.58 

11 

0.013 

0.026 

22 

86.99 

4 

0.005 

0.012 

23 

92.75 

1 

0.002 

0.004 


TOTAL 

10441 

1.000 

1.000 


D(10) - 17.81 pm 

D(20) - 23.04 nm 

D(21) = 29.73 nm 

D(30) - 27.69 pm 

D(31) - 34.48 /im 

D(32) - 40.01 pm 


1 Reproduced from specification sheet supplied by the manufacturer. 

2 Diameters traceable to NBS Part. #52577, 
accurate to ± 2 pm (± 3% for D > 70 pm) 
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(MAGL [10]) is calculated by thresholding the image at the optimum value (i.e., 90 T as specified 
by Ahlers), summing the pixel grey levels (GL) corresponding to specific particles as specified by 
segmentation, and dividihg by the total number of pixels per particle (Equation 2.5). 


MAGL = 


E iti GL{i,j) 

Pixel(t,j) 


(2.5) 


The PBG is determined by thresholding the image twice, once at 90 T, and the second, just below 
background (T*,). Referring to Fig. 2.12, the double threshold specifies the particle boundary 
gradient by: 

PBG = D d - D b , (2.6) 


where D& is the particle diameter at Tj,. With the above parameters, focus was specified for a 
volume centered on the focal plane of the transfer lens. First, a relation, constant with respect to 
focal volume, was determined for the MAGL with dependence on particle diameter, and second, 
the PBG was specified as a constant over the range of particle diameters specified by the MAGL 
criteria. 

In conclusion, Wiles developed a focus criteria for the LI/VPS. In his follow-up tests, the 
criteria defined a depth of focus which remained fairly constant when tested with the reticle and 
the polystyrene spheres (i.e., 52.5 fim as specified earlier). The prescribed depth of focus was 
approximately 400 microns. It should be noted, Wiles’ focus classification was determined and 
tested with the laser pulsing at 60 Hz. Thus, the focus criteria specified a depth of focus and 
classified particles based on grey level intensity from these operating conditions. 


2.2.5 Modifications 

The final goal of this research was the implementation of a particle sizing system capable of 
performing analysis on two- phase flow (e.g., aerosol sprays). The LI/VPS has been developed 
in stages; (1) Ahlers’ initial work, hardware and software setup, (2) Wiles’ work on system focus 
classification, and (3) the the current adaptation of the system to process truly dynamic particles 
in a real spray. To clarify the above statement, previous work by Ahlers and Wiles was performed 
with the LI/VPS operating in the continuous pulse mode (CPM), as opposed to the current work 
in the single pulse mode (SPM) (i.e., CPM suggests the imaging laser is pulsing at 60 Hz. in sync 
with the camera, and SPM implies the imaging laser is off until the video processor requires a new 
frame to process at which time the imaging laser is pulsed). The following discussion covers the 
reasoning and implementation of the SPM, and the adaptation of the previous work to function in 
the SPM. 

All previous work on the LI/VPS was done in the CPM, therefore the system had to be 
converted to the SPM. The reasoning for the conversion is shown in Fig. 2.13. The two graphs 
were taken with the system in the CPM; the only difference being the bottom particle is dynamic 
whereas the top particle is stationary. As shown, there is a significant reduction in intensity for the 
dynamic particle as opposed to the stationary particle. The above behavior is due to the camera 
tube’s ability to refresh between successive frames. In the CPM, the dynamic particle being frozen 
by the 10 ns laser pulse is present in the field of view for less than 16.67 ms (i.e., the time necessary 
to complete one field), but the static particle in the CPM shows greater intensity because of the 
cumulative effect of the particle blanking out the same area on the camera tube. The behavior 
being time-dependent implies the camera tube reaches a constant intensity after a sufficient amount 
of time. Because the software was developed for the system operating in the CPM, and all previous 
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a. - Static 92.75 pm Particle in the CPM. 


b. - Dynamic 92.75 pm Particle in the CPM 
Figure 2.13: Static vs. Dynamic Particle Representation in CPM 
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work was performed on static particles (i.e., particles which have motion but appear static to the 
system), the system had to be adapted to size dynamic particles. Revision to the system could be 
achieved by either changing the system software, or changing the system hardware. Figure 2.14 
shows the, MAGL vs. particle size, focus classification curves. As is shown, the ‘ dynamic y curve 
is less distinct than the 'static' curve. Because of the added ambiguities in the ‘ dynamic ' curve, a 
method had to be determined to simulate the behavior of the stationary particles for the dynamic 
particles. 

Because of the amount of work put into the development of the system software and the success 
of the focus criteria, a hardware modification was selected to accomplish the intensity contrast in 
dynamic particles. The SPM was found to exhibit the same characteristic intensity in the dynamic 
particles as found in static particles, in fact, the contrast between particle and background was 
greater. The SPM was accomplished by; (1) sending a trigger signal from the control computer 
to the LSC, (2) the LSC triggers the Na laser, (3) the laser pulses, and (4) the image processor 
grabs the frame just illuminated. The above procedure was accomplished by the development of 
a triggering circuit (APPENDIX B). The above procedure is then followed by normal program 
execution. The flow diagram in Fig. 2.15 shows the SPM integrated into the PSP with software 
modification. 

The software had to adapted to handle the SPM. As stated previously, the use of the SPM 
produced even greater contrast between the particle image and background. Because of the greater 
contrast, it was necessary to redetermine the focus criteria. Using the procedure outlined by Wiles 
[10] (Section 2.2.4), the MAGL curve and the PBG criteria were determined in the SPM. MAGL 
curves for both the CPM and the SPM are represented in Fig. 2.16. As shown in the figure, the 
larger particles show greater contrast whereas the smaller particles contrast is unaffected by the 
SPM. The focus criteria was determined for both the 500X and 1000X lens. The LI/VPS, at this 
point, was capable of performing size measurements in a two-phase flow. 

2.2.6 Software Updates 

With PSP performing analysis on two-phase flows, the software had to be updated to allow for 
varying conditions in the measurement analysis. Parameters, such as the sizing window specifi- 
cations, output destination, etc,, were queried for before processing each time the program was 
executed and others, such as lens magnification, were set by changing the FORTRAN code. A 
menu type of setup (Fig. 2.17) was adopted to minimize setup time and to aid the operator in 
determining the most appropriate sizing conditions (APPENDIX C.l). 

In aerosol 6prays, the mean diameters (APPENDIX D) determined from the count vs. drop-size 
data are the most common method of characterization. Characterization by mean diameters is mis- 
leading when a single mode (i.e., Gaussian distribution) is not the case, therefore the actual count 
vs. drop-size distribution is also used to characterize aerosol sprays. Because of the aforementioned 
reasoning and the unavailability of a suitable graphics package for the LI/VPS, a graphic algorithm 
was developed. The algorithm was coded into a FORTRAN subroutine (APPENDIX C.2 ) for the 
PSP with a DEC VT240 terminal for graphic simulation (Fig. 2.18(a)) and a DEC LA75 printer 
for hard-copies (Fig. 2.18(b)). 

2.3 Spray Test Facility 

Figure 2.19 shows the configuration of equipment for the spray characterization tests. The tests 
were performed in the horizontal direction due to the positioning of the sizing instrumentation. 
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AAA Points Used to Determine MAGL for Dynamic Particles 

MAGL Parameter for Dynamic Particles 

MAGL Parameter for Stationary Particles 


Diameter (pm) 


Figure 2.14: Comparison of the MAGL Parameter for Dynamic and Stationary Particles in CPM 






Figure 2.15: Flow Diagram for the Particle Sizing Program in SPM 
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SETUP 

PARTICLE SIZING PROGRAH (ver. 4) 


- - PROCESSING 

(A) DYNAHC Type of Processing 

(B) TES Focus Criteria 

(C) AUTO Type of Fraee Advance 

(0) PARTCL Processing Limit 

(E) ( 1000) Limiting Value 

(F) REJECT Boundary Particles 

OUTPUT 


( 0 ) 
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(J) 

(N) 

(O) 
(S) 
(U) 


OPTIONS 

(STATIC/DYNAHIC) 

(YES/NO) 

(AUTO/SINGLE) 

(TIME/FRAME/PARTICLE) 

( seconds/ frames/par tides) 
(PROCESS/REJECT) 

OPTIONS 


YES 

NO 

NO 

YES 


General Results (to PRINTER) 
WRITE TO FILE (YES/NO) 
Average Particle sits data -■ 

Group Breakdown data 

Per Frame data : 


(YES/NO) 

(K) FILE HEADER (4 lines) 
FXLEi (TEHP01 ) . OUT 


(L) 
•/ 

(M) 


Group Start •( 5) (0) 

X Window Start - ( SO) 

Y Window Start ■ ( SO) 

Threshold * ( 90) (V) 


GENERAL 

Group Width 


Lens 


HIGH 


OPTIONS 
•( 5.0) 

(R) 

(T) 

<*») 


FILE! (TEMP01 ) . DAT 

(P) W of Groups • ( 68) 

X Window Width - (450) 

Y Window Width - (4S0) 
Markers - NO (YES/NO) 


(X) to exit SETUP menu or (Z) to begin Partlcal Siting Program .... 
Enter Letter to change specific Parameter 77 


Figure 2.17: PSP Setup Page 
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DEC VT240 Terminal for Screen Emulation. 



b. - DEC LA75 Printer for Hard-copys. 
Figure 2.18: PSP Graphic Package 
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Figure 2.19: Equipment Schematic for Instrument Comparison 
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The experimental apparatus was situated on a Newport Research optical table equipped for isola- 
tion. The building ventilation system was used to draw off the aerosol spray after analysis. The 
spray characterization tests were performed on an air-assist nozzle. 

2.3.1 MOD-1 Nozzle 

Figure 2.20 shows the MOD-1 nozzle as supplied by NASA Lewis Research Center. The nozzle 
is of the atomizer type and a prototype of the nozzle proposed to be used in the NASA Altitude 
Wind Tunnel to simulate various cloud structures in icing studies. Variation of the drop-size in 
the aerosol spray produced by the nozzle is obtained by varying the input air and water pressures. 
The water is introduced into a 1.81 inch long by 0.368 inch diameter mixing chamber through a 
0.0155 inch orifice. The air is introduced into the outer wall of the mixing chamber through twelve 
0.125 inch holes. After mixing, the aerosol is expelled from the mixing chamber through a 0.125 

inch orifice. 

2.3.2 Air and Water Supply System (AWSS) 

As shown in Fig. 2.21, the AWSS was constructed to supply air and water to the MOD-1 nozzle 
with the exception of the LI/VPS optics purge. The air for the AWSS is supplied by twin 100 
hp Ingersoll-Rand turbine compressors with a delivery rate of 800 SCFM at 120 psig. Because of 
the high water pressure necessary for the MOD-1 nozzle, a Brunswick 20.5 liter pressure vessel 
was filled with water and pressurized by the supply air or for higher pressures by a regulated high 
pressure N 2 bottle. After pressurization, the water was filtered by an ADKIN spool filter. The 
nozzle air and water supply was regulated by a WATTS Model 2235 pressure regulator and a Cole- 
Parmer Model PR004-FM 044-4 0G flowmeter, respectively. Connection lines in the supply system 
were YELLOW JACKET Model WPP0031A charging hose (500 max. psi.). The LI/VPS optics 
purge used a regulated high pressure Nj bottle for a constant positive flow from the lens cover to 
avoid contamination. 

2.3.3 Water Flowmeter Calibration 

The Cole-Palmer flowmeter was factory calibrated. The calibration was verified by collecting and 
weighing the water passing through the flowmeter. The water was weighed on a HOWE model 
#3074131 balance scale. Twelve flow rates were measured with three samples collected at each 
flow rate. The experimental data and factory calibration data are presented in Table E13.1 with 
graphical representation shown in Fig. E.13.1 (APPENDIX E). 


2.4 Digital Pressure Acquisition 

The digital pressure system (DPS) was developed to monitor the essential input conditions of the 
MOD-1 nozzle. The DPS consists of two OMEGA Model PX304-150AV pressure transducers, a 
DEC AXV11- C analog to digital (A/D) converter board, the PDP- 11/73 micro- computer hosting 
the above A/D board, and a PDP RT-U software package written to access the A/D board and 
store or display the resulting pressures. 
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b. - MOD-1 Components. 
Figure 2.20: MOD-1 Nozzle 
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Figure 2.21: Air and Water Supply Schematic 
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2.4.1 Pressure Transducers 


The OMEGA pressure transducers (Fig. 2.22) are bridge type strain gage transducers. The bridge 
excitation voltage was 10 VDC supplied by a Hewlett-Packard (Model Harrison 6200B) d.c. power 
supply with a bridge output of 0 to 100 mVDC. The transducers are specified to have an operating 
range of 0 to 150 psia with ± 0.75 psi accuracy. 

2.4.2 A/D converter board 

The DEC AXV11-C analog- to-digital converter board was installed in the back-plane of the PDP- 
11/73 microcomputer. The AXV11-C board has 12 bit digital resolution, supports up to 16 single 
analog input signals or 8 differential signals, A/D conversion by program, external clock, or real- 
time clock, and 1, 2, 4, and 8 (i.e, 10, 5, 2.5, and 1.25 volts) programmable gain settings. As 
recommended by the manufacturer, the 8 channel differential option was chosen to maximize analog 
to digital conversion, due to the 100 mV range supplied by the pressure transducers. 

2.4.3 Analog- to- Digital Conversion 

The transducer voltage signal is converted to a digital value available to the LI/VPS operator. An 
interface box (Fig. 2.23) was constructed to utilize the full capabilities of the AXVll-C board. The 
interface box has 8 A/D input ports and 2 D/A output ports using BNC connectors. The interface 
box is linked to the AXVll-C board by R.S232 cable and connectors. The pressure measurements 
are made available to the analyst through the PDP-11/73 microcomputer. The RT-11 software 
package, w’ritten in FORTRAN subroutine form (APPENDIX C.3), allows for real-time pressure 
monitoring with storage and averaging capabilites for the duration of the main calling program. 
The A/D converter is programed for a gain setting of 8 (i.e., an effective analog input range of 0 
to 1.25 volts) to optimize A/D conversion of the pressure tramsducer output range of 0 to 100 mV. 

2.4.4 Digital Pressure System Calibration 

The pressure transducers were calibrated for various static pressures by pressurizing the transducers 
and reading the A/D output after a steady equilibrium state had been attained. A laboratory grade 
test gage was used to measure the “ standard" pressure. The test gage, with a range of 0 to 160 
psig, was calibrated using an American Steam Gage Co. deadweight pressure gage tester. With 
the pressure tramsducer ’s specified input pressure range of 0 to 150 psia, the calibration data was 
taken within a range of 0 to 110 psig ( 14.05 to 124.05 psia). The atmospheric pressure at the time 
of the calibration run was measured to be 727.29 mm Hg. or 14.05 psia from a Precision Thermo & 
Inst Co. model #Z769 barometer. The experimental data is presented in Tables F14.1 and F14.2 
with graphical representation shown in Figs. F14.1 and F14.2 (APPENDIX F). 

2.5 Experimental Procedure 

With system performance and verification as the basis for comparison, equivalent sampling was 
required. As discussed earlier, the P/DPA and the LI/VPS use different methods of particle sizing 
(i.e., temporal vs. spatial), but each instrument uses a probe volume for data collection. Therefore, 
system comparison was dependent on spray density, droplet size range, and user designation of the 
measurement volumes (i.e., the P/DPA ’s crossed-beam intersection volume, specified by 
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Model # - PX 304-150A V 


SPECIFICATIONS 
Excitation: 10 VDC 
Output: 0 to 100 mV 
Sensitivity: 10 mV/V ±1% 

Input Impedance: 1200 ohm 
Output Impedance: 500 ohm 


PERFORMANCE 
Accuracy: ±0.5% full scale 
Zero Balance: ±2.0% full scale 
Operable Temperature Range: 
-29 to 60° C 


a. - OMEGA Pressure Transducer Data Sheet. 



b. - OMEGA Pressure Transducers. 
Figure 2.22: OMEGA Pressure Transducers 
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Figure 2.23: A/D Connector Box Schematic 
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the transmitter lens chosen and beam diameter, vs. the LI/VPS focus volume, specified by the 
imaging optics and software). 

The procedure for overlapping the probe volumes is described in reference [14]. Figure 2.24 
is included to show the scattered light, from drops generated by the VOAG passing through the 
crossed-beam intersection volume, as seen by the LI/VPS. 


2.5.1 Verification Tests 

The P/DPA and the LI/VPS probe volumes for the verification tests were specified as follows; the 
P/DPA transmitter lens with the 495 mm focal length and 25 mm beam separation formed a probe 
volume with an approximate 160 fim waist diameter, and for the LI/VPS, the 1000X lens specifies 
a 400x400x140 /xm 3 volume with software selectable field of view for a 160x160x140 /im volume 
(Fig. 2.25). 

With the above configuration, the P/DPA and the LI/VPS were tested using a TSI Model 3450 
Vibrating Orifice Aerosol Generator (VOAG). Operating conditions of the VOAG were varied to 
generate a size range of particles, 19.8 to 99.6 /xm (Table 2.2). 


TEST 

(#) 

Table 2.2: Verification Test Conditions 
ORIFICE DISTURBANCE WATER 

DIAMETER FREQUENCY FEED RATE 
(ixml (Hz.) (cm 3 /min) 

THEORETICAL 

DIAMETER 

(/xm) 

v ” J- 

1 

10 

330.4 

0.080 

19.8 

2 

20 

100.2 

0.139 

35.5 

3 

20 

79.2 

0.139 

39.0 

4 

20 

62.5 

0.139 

41.5 

5 

20 

51.6 

0.139 

44.2 

6 

20 

41.6 

0.139 

47.5 

7 

50 

30.1 

0.590 

85.6 

8 

50 

25.5 

0.590 

90.4 

9 

50 

19.0 

0.590 

99.6 


Each particle size generated either in single stream form or using the dispersion cup (Fig. 2.26)to 
generated a spray was measured using the P/DPA and the LI/VIPS system. The TSI droplet 
diameter (Dj) was calculated using the TSI theoretical equation (2.7); 

MSP <27) 

where q is the liquid flow rate and f is the disturbance frequency. Results of the tests are presented 
in Section 3.2. 
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Dispersed Droplet 



a. - VOAG Dispersion Cup. 



b. - TSI Vibrating Orifice Aerosol Generator. 
Figure 2.26: Verification Test Apparatus 
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2.5.2 Spray Comparison 

With the spray density and particle size range depending on the nozzle conditions, the benchmark 
tests were performed for two specific cases. Inlet nozzle conditions are shown in Table 2.3. 

Table 2.3: Comparison Test Conditions 

CASE I CASE II 

Pressure (water) 115 psia 105 psia 
Pressure (air) 45 psia 65 psia 


For each case, a sample was taken on centerline two feet downstream from the nozzle with suc- 
ceeding samples taken radially in 0.5 inch increments to the outer edge of the spray plume. 

To avoid undue comparative data reduction, the P/DPA and LI/VPS were matched in approxi- 
mate probe volume size, as previously stated, and appropriate particle size range. Assuming nozzle 
conditions were steady state, preliminary setup of the P/DPA and the LI/ VPS was performed to 
optimize instrument operation. The results of the analysis are presented in Section 3.3. 
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Section 3 


PRESENTATION AND DISCUSSION OF RESULTS 


This section will present the results of the LI/VPS calibration tests including a comparison with 
previous calibration tests, the verification tests with the VOAG, and the comparison tests using 
the MOD-1 nozzle. The major concern of these results is the accuracy of the sizing measurements 
with secondary interest in the comparability of the LI/VPS and the P/DPA. 


3.1 LI/ VPS Calibration Results 

As was stated previously, the LI/VPS had to be recalibrated due to the replacement of the vidicon 
camera tube. With the new vidicon tube, the MCF became approximately 2.1 /xm/pixel (i.e., 
for the 500X lens), as opposed to the previous factor of 1.8 /xm/pixel [7,10], for the old camera 
tube. The new vidicon tube, therefore, reduced the LI/VPS measurement resolution. The above 
is mentioned to explain the increased error in determining the smaller particle sizes for the 500X 
lens, as well as the reasoning for the calibration of the 1000X lens. The following calibration results 
specify the MCFs for the 500X and the 1000X lens. Results of previous calibration tests using the 
calibration reticle have been compared to the new calibrations. 

Using the procedure described in Section 2.2.3, the Equations (3.1) thru (3.4) represent the 
MCFs as functions of x and y location for the two lens; 


the xMCF for the 500X lens; 


- 04 for 50 < x < 100 


MCF(y) = 


' 2.21 + y 
2.20 + y 
2.16+ y 
2.16 + y 
2.16 -y 
2.11 -y 
2.10 — y 
. 2.07 - y 


* 0.803F 

* 0.290F 

* 0.679F 

* 0.442E 

* 0.947F 

* 0.306F 

* 0.124F 

* 0.1351? 


o 

1 

for 

100 

i 

o 

for 

150 

r- 

o 

1 

for 

200 

o 

1 

for 

250 

1 

o 

for 

300 

-03 

for 

350 

- 03 

for 

400 


< x < 150 

< x < 200 

< x < 250 

< x < 300 

< x < 350 

< x < 400 

< x < 450, 


(3.1) 


39 



the yMCF for the 500X lens; 


< 


MCF(y) = { 


2.10 -y* 0.1832? - 03 for 50 < x < 100 

2.11 -y* 0.240£ - 03 for 100 < x < 150 

2.12 - y * 0.314F - 03 for 150 < x < 200 

2.13 - y * 0.313F - 03 for 200 < x < 250 

2.15 - y * 0.3972? - 03 for 250 < x < 300 

2.18 - y * 0.4842? - 03 for 300 < x < 350 

2.19 - y * 0.5052? - 03 for 350 < x < 400 

2.18 — y* 0.509.E - 03 for 400 < x < 450, 


the xMCF for the 1000X lens; 


/ 


MCF(y) = { 


0.977 + y * 8.09 E - 05 for 50 < x < 100 
0.974 + y * 2.60 E - 05 for 100 < x < 150 
0.967 - y * 8.12 E- 07 for 150 < x < 200 
0.961 + y * 4.732? - 06 for 200 < x < 250 
0.961 - y * 5.46 E - 05 for 250 < x < 300 
0.948 -y* 3.72 E - 05 for 300 < x < 350 
0.943 - y * 6.80 E - 05 for 350 < x < 400 
0.920 - y * 2.58 E - 05 for 400 < x < 450, 


and the yMCF for the 1000X lens; 


MCF(y) { 


0.977 - 
0.981 - 
0.981 - 
0.990 - 
1.000 - 
1.014- 
1.027 - 
1.029- 


y * 9.17 E - 05 for 50 < x < 100 
y * 1 .242E7 - 04 for 100 < x < 150 
y * 1.192? - 04 for 150 < x < 200 
y * 1.632? — 04 for 200 < x < 250 
y * 1.96 E — 04 for 250 < x < 300 
y * 2.19 E - 04 for 300 < x < 350 
y * 2.63 E - 04 for 350 < x < 400 


y * 2.69 E - 04 for 400 < x < 450. 


(3.2) 


(3.3) 


(3.4) 


With the above equations, a software algorithm was setup in subroutine form to determine the cor- 
rection factors as functions of particle location and for the magnification lens installed (APPENDIX 

c -4)- 

Figures 3.1 - 3.4 show the variation of the MCFs with respect to x and y location. The similarity 
in Figs. 3.1 and 3.3, as well as the similarity in Figs. 3.2 and 3.4 show the MCFs’ variation is mainly 
due to the geometric non-lmearities in the vidicon tube. The procedure developed to determine 
the MCFs as functions of both x and y screen location is easy to use, straight-forward, and not 
time consuming. The implementation of the MCFs in PSP is easily facilitated by the use of the 
FORTRAN subroutine format. 

The following comparison represents LI/VPS accuracy studies by this investigator and the 
previous investigators [7,10]. The basis for the comparison was the utilization of the calibration 
reticle with the 500X lens. Table 3.1 shows the results for the 500X lens by this investigator. Table 
3.2 represents the equivalent results for the 1000X lens under similar test conditions. 
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Figure 3.1: Magnification Correction Factor Behavior 
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Figure 3.3: Magnification Correction Factor Behavior 
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Table 3.1: Calibration Accuracy Test 
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Table 3.2: Calibration Accuracy Test 


Table 3.3 shows the average percent error for the above calibration accuracy tests with the previous 
work of Ahlers [7] and Wiles [10]. A comparison of the average % error for the three accuracy tests 
performed on the 500X lens shows a decrease in the % error from the one- dimensional MCF test 
(i.e., 4.04% error) to the two-dimensional MCF tests (i.e., for Wiles - 2.73% error and for this work 
- 3.96% error). The % error values for the test performed on the 1000X lens show an increase in 
LI/VPS accuracy for all the particles measured by the 500X lens tests. The inclusion of the 5.29 
H m particle in the analysis shows an increased sizing range, as opposed to previous tests. 

The following results represent the initial method used to compare the P/DPA and the LI/VPS. 
As specified earlier, the probe volumes of the two instruments were overlapped, and due to the 
steady state operation of the VOAG, samples by both instruments were assumed to be nearly 
identical. Two separate cases were performed to verify instrument operation and accuracy. The 
first case was performed with the VOAG producing a steady single stream of drops which passed 
through the concurrent probe volumes, and secondly, the dispersion cup (Fig. 2.26) was utilized 
to produce a spray of monodisperse droplets which randomly pass through the concurrent probe 
volumes. Nine separate tests were performed for each case with the instrument results represented 
in Figs. 3.5 thru 3.13 for the case without the dispersion cup, and Figs. 3.14 thru 3.22 for the 
case with the dispersion cup. Figures 3.23 and 3.24 show the TSI theoretical diameter, and the 
arithmetic mean diameters from the LI/VPS and the P/DPA distributions as functions of test 
number. Data in Table 3.4 has been plotted in Fig. 3.23 and 3.24 with the standard deviation 
(SD) also shown. The arithmetic mean diameters of the LI/VPS and the P/DPA agree, on the 
most part, with each other and the theoretical expected diameter within ± 2.6 /xm. The SD of 
the samples is shown to illustrate the monodisperse behavior of the VOAG and the ability of the 
LI/VPS and the P/DPA to measure the monodisperse aerosol spray. The highest SD (i.e., 1.109 
/im) determined for the LI/VPS is shown in CASE II - Test 5, and for the P/DPA, the highest SD 
(i.e., 2.073 /xm) is shown in CASE I - Test 1. 

Referring to Table 3.4, the first test in both cases show the maximum SD for P/DPA. The 
arithmetic mean diameters, 20.5 /xm for CASE I and 21.5 /xm for CASE II, are within 2.0 /xm of 
the expected diameter, 19.8 /xm. The SD of the samples may be higher than the rest, due to the 
high density of drops passing through the P/DPA probe volume. This phenomena was especially 
noticeable in CASE II test runs where the dispersion cup was used. As was expected, the SD for 
most of the tests increased from CASE I to CASE II. The above behavior was expected, due to 
the increase in number of drops passing through the edges of the probe volumes. 

3.2 Results For the MOD-1 Nozzle Comparison 

The following results represent a comparison of the LI/VPS and the P/DPA in side-by-side bench- 
mark tests performed on a NASA MOD-1 atomizing nozzle. As previously stated, two cases (i.e., 
variation in the operating conditions of the nozzle) were studied. For each case, eight data runs (i.e., 
a data run was performed on the centerline, two feet down-stream from the nozzle with succeeding 
data runs performed at one-half inch increments radially outward to the edge of the dispersion) 
were performed by the LI / VPS and the P/DPA using a procedure similar to the VOAG analysis. 
Figures 3.25 - 3.32 and Figs. 3.33 - 3-40 are the results from the P/DPA and the LI/VPS for CASE 
I (i.e., nozzle conditions: Air pressure = 65 psia and Water pressure = 105 psia.) and CASE II 
(i.e., nozzle conditions: Air pressure = 45 psia and Water pressure =115 psia), respectively. 
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Table 3.2: Calibration Accuracy Test 


CALIBRATION RETICLE : RR-50-3.0-0.08- 102-CF - *114 


For the 500X Lens. For the 1000X Lens 


PART. 

(#) 

DIAMETER 

(pm) 1 

Ahlers’ 
TEST (7] 
Avg. % Error 

Wiles’ 
TEST [10] 
Avg. % Error 

Current 

TEST 

Avg. % Error 

Current 

TEST 

Avg. % Error 

1 

5.29 

0.00 

0.00 

0.00 

33.65 

2 

6.81 

17.75 

16.06 

26.58 

17.77 

3 

8.98 

7.66 

7.40 

8.69 

4.23 

4 

11.93 

12.12 

7.53 

2.77 

5.28 

5 

17.20 

3.64 

1.52 

10.41 

4.59 

6 

21.33 

3.33 

2.04 

5.77 

3.89 

7 

23.90 

2.05 

2.19 

5.61 

3.10 

8 

26.71 

3 01 

3.89 

3.03 

4.16 

9 

31.11 

5.06 

1.66 

2.09 

3.28 

10 

34.17 

3.27 

1.66 

3.43 

2.25 

11 

37.07 

2.34 

1.21 

1.79 

1.54 

12 

40.47 

1.04 

1.91 

1.66 

4.13 

13 

42.71 

1.50 

2.50 

1.19 

3.39 

14 

47.37 

3.96 

1.52 

1.41 

0.78 

15 

50.39 

3.44 

0.94 

0.77 

1.51 

16 

52.50 

2.06 

1.21 

1.03 

0.65 

17 

56.23 

1.35 

1.09 

0.44 

2.22 

18 

60.70 

0.91 

1.33 

1.04 

3.23 

19 

67.04 

4.50 

0.83 

0.66 

0.60 

20 

73.48 

3.02 

1.25 

3.84 

0.46 

21 

80.58 

2.50 

0.76 

4.28 

0.53 

22 

86.99 

1.45 

0.88 

1.80 

1.41 

23 

92.75 

0.92 

0.67 

0.70 

0.38 

AVERAGES: 

4.04 

2.73 

3.96 

4.48 


1 Diameters traceable to NBS Part. *52577, 
accurate to ± 2 /un (± 3% for D > 70 pm) 
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Figure 3.5: VOAG Verification w/o Dispersion Cup 
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Vibration Frequency = 100.2 kHz 

Figure 3.6: VOAG Verification w/o Dispersion Cup 
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Figure 3.7: VOAG Verification w/o Dispersion Cup 
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Figure 3.8: VOAG Verification w/o Dispersion Cup 
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Diameter of Orifice = 20 //m 
Liquid Feedrate = 0.139 cm 3 /min. 
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Figure 3.9: VOAG Verification w/o Dispersion Cup 
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Liquid Feedrate = 0.139 cm 3 /min. 
Vibration Frequency = 41.6 kHz 

Figure 3.10: VOAG Verification w/o Dispersion Cup 
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Figure 3.11: VOAG Verification w/o Dispersion Cup 
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Figure 3.12: VOAG Verification w/o Dispersion Cup 
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Figure 3.13: VOAG Verification w/o Dispersion Cup 
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Theoretical Diameter = 19.8 /xm 
Diameter of Orifice = 10 ^m 
Liquid Feedrate = 0.08 cm 3 /min. 
Vibration Frequency = 330.4 kHz 

Figure 3.14: VOAG Verification w/ Dispersion Cup 
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Figure 3.15: VOAG Verification w/ Dispersion Cup 
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Figure 3.16: VOAG Verification w/ Dispersion Cup 
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Figure 3.18: VOAG Verification w/ Dispersion Cup 
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Figure 3.19: VOAG Verification w/ Dispersion Cup 
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Figure 3.20: VOAG Verification w / Dispersion Cup 
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Figure 3.21: VOAG Verification w / Dispersion Cup 
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70 








f» n 



658 


20.6 

Dianeter (mcroneters) 


Most ProWble Du: 7.1 
An thnetic Mean (D 10 ): . I 

Am Hew <D 20 = 13.6 
Volume Hew (D 30 ): 14.9 
Sauter Hew ( 032 ): 18.2 


Corrected Count: 12290 


File: HAS 4 22 .MT 
Date: 08 - 0^-1987 


Atnp: 12258 
Tine: 10:01:11 


Total Count: 
Sun Tine: 


9375 

5.148 


«< 

seconds 



371 


Velocity Hew 
RMS velocity 


Velocity 


11.2 

(neters/second) 


20.0 


PhS= 

Ovjc 

Und= 

Hax= 


0 

270 
1109 
1 


Sat: 

•ia= 

Ovf= 

VeU 


14.71 

2.8 


Total Bad = 2100 


4 

590 

0 

122 


boLdtv 


1109 


a, P/DPA Results 


w 

*- 

z 


50 


Spatial Distribution 


37 


25 


12 


Distribution Mode Dio* 
Ari thnetic Hean Dio. 
Surface Hean Dia. 
Volune (Hass) Hean Dia. 
Sauter Hean Dia. 


Total 


X 

7.0 

(D10) = 

9.8 - 

(D20) = 

11.2 

(D30) = 

12.8 

(D32) = 

16.5 

Count * 

500 — 



10.8 20.5 

DIAMETER (nicrons) 


30.3 


40.0 


b. LI/VPS Results 
(CASE I) 

Test Conditions: Radial Position = ^ in. 
Air Pressure = 65 psia 
Water Pressure = 105 psia 
Water Flowrate = 0.038 gal/min. 
Axial Position from Nozzle = 2 ft. 


Figure 3.26: MOD-1 Nozzle Comparison 


71 




r» o 



693 


28.6 

Di am ter (micrometers) 


Host Probable Dia: 7.1 
Arithmetic Mean (DIB): 11,9 
Area Mean (D29): 13.4 
Volume Mean (038): 14.9 
Sauter Mean (D32): 18.4 


Corrected Count: 12282 


File: NAS4J3.DAT 
Date: 88-84-1987 


Atmp: 11971 
Time: 18:84:56 


Total Count: 
Run Tine: 


18815 

4.784 


«< 

seconds 



351 


Velocity Neon 
RMS velocity 


13.7 

3.1 


2.4 11.2 

Velocity (meters/second) 


28.8 



605 


a. P/DPA Results 


a 

h 

Z 


60 


45 


30 


15 


1 r 

Spatial Distribution 


Distribution Mode Dia. 
Arithmetic Mean Dia. 
Surface Mean Dia. 
Volume (Hass) Mean Dia. 
Sauter Hean Dia. 


a 

5.9 

(DiO) « 

9.5 — 

(D20) 3 

10.9 

(D30) 3 

12.4 

(DS) 3 

16.2 

Count 3 

506 — 


1.0 



10.8 


20.5 


30.3 


40.0 


DIAMETER (microns) 
b. LI/VPS Results 
(CASE I) 

Test Conditions: Radial Position = 1 in. 
Air Pressure = 65 psia 
Water Pressure = 105 psia 
Water Flowrate = 0.038 gal/min. 
Axial Position from Nozzle = 2 ft. 


Figure 3.27: MOD-1 Nozzle Comparison 


72 





830 


20.6 

Dianeter (nicroneters) 


File! NAS4 24. DAT 
Date: 08-04-1987 


Host Pro table Di a- 7.1 
Anthnetic dean (D10): 10.8 
Area dean (D20)= 12.3 
Volune Mean (D30>= 13.7 
Sauter dean C D32 ) = 17.2 


Corrected Count: 13483 


A tap: 12117 
Tim: 10:08:38 


Total Count: 1 
Run Tine: 


3.236 


«< 

seconds 



Velocity 


11.2 

(neters/second) 


Velocity dean = 12.48 
RMS velocity = 3.33 


Phs= 0 $at= # 0 

Ovr= 282 Dia= 647 

Und= 577 0 

dax= 1 Vel= 34 

Total Bad = 1541 


EC 

BOBKDTV 


647 


a. P/DPA Results 


Spatial Distribution 


Distribution Mode Die* 
Arithnetic Mean Dia. 
Surface Mean Dia. 
Vo 1 une (Hass) Mean Dia. 
Sauter Mean Dia. 



> 

7.0 

(D10) = 

10. i — 

(D20) = 

11.3 

(D30) = 

12.8 

(D32) = 

16.2 

Count - 

503 — 


I.ltll.I l-l_ i 

30.3 


40.0 


DIAMETER (nicrur.s) 


b. LI/VPS Results 
(CASE I) 

Test Conditions: Radial Position = 1^ in. 
Air Pressure = 65 psia 
Water Pressure = 105 psia 
Water Flowrate = 0.038 gal/min. 
Axial Position from Nozzle = 2 ft. 

Figure 3.28: MOD-1 Nozzle Comparison 


73 





illllllliLl..ii.i..liii.. L... 

10.8 20.5 30.3 

DIAMETER (*icront) 
b. LI/VPS Results 
(CASE I) 

Test Conditions: Radial Position = 2 in. 
Air Pressure = 65 psia 
Water Pressure = 105 psia 
Water Flowrate = 0.038 gal/min. 

Axial Position from Nozzle = 2 ft. 


Figure 3.29: MOD-1 Nozzle Comparison 


74 












DIAMETER (xicrons) 
b. LI/VPS Results 


(CASE I) 

Test Conditions: Radial Position = 2| in. 
Air Pressure = 65 psia 
Water Pressure = 105 psia 
Water Flowrate = 0.038 gal/min. 
Axial Position from Nozzle = 2 ft. 

Figure 3.30: MOD-1 Nozzle Comparison 


75 





DIAMETER (nicrons) 


b. LI/VPS Results 


(CASE I) 

Test Condtions: Radial Position = 3 in. 
Air Pressure = 65 psia 
Water Pressure 105 psia 
Water Flowrate = 0.038 gal/min. 
Axial Position from Nozzle = 2 ft. 


Figure 3.31: MOD-1 Nozzle Comparison 


76 




a. P/DPA Results 



1.0 9.5 18.0 26.5 


DIAMETER (microns) 
b. LI/VPS Results 
(CASE I) 

Test Condtions: Radial Position = 3^ in. 
Air Pressure = 65 psia 
Water Pressure = 105 psia 
Water Flowrate = 0.038 gal/min. 
Axial Position from Nozzle = 2 ft. 

Figure 3.32: MOD-1 Nozzle Comparison 


77 


is 




File: HAS4JI.MT 
Date: 88-94-1987 


(I tup: 12164 
Tim: 99:91:25 


Host Probable Du: 19.9 
Arithmetic Hean (D19): 19 
Area Mean (Dll): 23 
Volume Hean < D39)= 27.4 

132): “ * 


Sauter Hean ( 3 


33.2 


Corrected Count: 14374 


Total 


Count: 19991 
i Tine: 2.1 


999 


«< 

seconds 



a. P/DPA Results 


75 


55 


37 


18 


0 +•* 

2.0 


Spatial Distribution 


Distribution Hodo Dia. 

* 

6.0 

Arithmetic Horn No. 

(DAO) * 

11.9 

Surface Hean No. 

(D20) 3 

15.3 

Volume (Hoss) Hean Dia. 

(D30) * 

19.6 

Sauter Hem Dio. 

(D32) 3 

32.5 

1 Total 

Count 3 

501 





21.5 41.0 

DIAMETER (microns) 


60.5 


80.0 


b. LI/VPS Results 
(CASE II) 

Test Conditions: Radial Position = Cl 
Air Pressure = 45 psia 
Water Pressure =115 psia 
Water Flowrate = 0.Q6 gal /min. 
Axial Position from Nozzle = 2 ft. 


Figure 3.33: MOD-1 Nozzle Comparison 


78 






a. P/DPA Results 



DIAMETER (nicrons) 
b. LI/VPS Results 
(CASE II) 

Test Conditions: Radial Position = ^ in. 
Air Pressure = 45 psia 
Water Pressure =115 psia 
Water Flowrate = 0.06 gal/min. 
Axial Position from Nozzle = 2 ft. 

Figure 3.34: MOD-1 Nozzle Comparison 


79 





Diameter (micrometers) 

File: NAS4 33.DAT ft tup: 12258 

late: 98-94-1987 line: 09:07:31 


92 

Most Probable Dia: 8.5 
Arithmetic Mean ( D10) = 16.2 
Area Mean ( D28) = 19.3 
Volume Mean (D30): 22.8 
Sauter Mean (D32): 31.7 


Corrected Count: 13958 

Total Count: 9504 «< 

Run Time: 2.226 seconds 



Velocity (meters/second) 


380 Velocity Mean : 16.74 
RMS velocity : 4.13 


Phs= 0 

0 vr= 301 

Und= 651 
Nax= 5 

Total Bad 


Sat= 
Dia: 1( 
Ovf: 
Vel= 

= 2561 


1600 


a. P/DPA Results 



2.0 21.5 41.0 60.5 80.0 


DIAMETER (microns) 

b. LI/VPS Results 
(CASE II) 

Test Conditions: Radial Position = 1 in. 
Air Pressure = 45 psia 
Water Pressure = 115 psia 
Water Flowrate = 0.06 gal/min. 
Axial Position from Nozzle = 2 ft. 

Figure 3.35: MOD-1 Nozzle Comparison 


80 





30.9 

Dianeter (mcroneters) 


9)7 


Host Probable Du- 9 
Arithmetic Hean (D10)e 15.4 
Area Mean (D20): 17.8 
Volume He an (D39>= 20.3 
Sauter He an < 032 ) = 26.7 


Corrected Count: 13066 


File: NAS4J4.DAT 
Date: 08-0^-1987 


Atnp: 12258 
Tine: 09:11:31 


Total Count: 
Run Tine: 


9113 «< 

3.572 seconds 



356 Velocity He an = 14.2 
IMS velocity = 4.2 


17.7 30.6 

Velocity (neters/second) 


Phs= 

9 

$at= 

Ovr= 

519 

»ia= 2 

Und= 

352 

0vf= 

Max? 

2 

Vel= 

Total Bad 

= 3918 


9 

2136 

4 

5 


QD 


2136 


a. P/DPA Results 



2.0 21.5 41.0 60.5 80.0 


DIAMETER (microns) 
b. LI/VPS Results 
(CASE II) 

Test Conditions: Radial Position = lj in. 
Air Pressure = 45 psia 
Water Pressure =115 psia 
Water Flowrate = 0.06 gal/min. 
Axial Position from Nozzle = 2 ft. 


Figure 3.36: MOD-1 Nozzle Comparison 


81 





8t>7 


30.9 

Dianeter (mcroweters) 


Host Probable Du* 81 
Arithmetic Hean vDl0): H.4 
Area Mean (D20): 16. 6 
Volume Hean (D30): 19.1 
Sauter Hean (D32): 25.1 

Corrected Count: 11378 


rile: NASA 35. DAT 
Date: 08-W-198? 


A tup: 12258 
Tine: 09:13:28 


Total Count: 
Run Tine: 


8023 «< 

5.259 seconds 



321 Velocity Hean = 11.78 
RHS velocity : 3.86 


4.8 17.7 

Velocity (meters/second) 


30.6 


Phs= 0 Sat: 0 

0w= 674 Dip 3169 

Und: 276 Ovf: 8 

Hax: 1 Del: 25 

Total Bad : 4153 BOUHSDTV 


!IZB 


3169 


a. P/DPA Results 


a 

t- 

z 

3 


65 


48 


32 


16 


1 r 

Spatial Distribution 


Distribution Node Dia. 
Arithmetic Hean Dia. (DIO) 
Surface Hean Dia. (D20) 
Volute (Hass) Hean Dia. (D30) 
Sauter Hean Dia. (D32) 


Ul 


III llllllillJll l ii. uJUI . 


Total Count = 


1.0 


15.8 


30.5 


45.3 


6.6 

11.9 

14.1 

16.3 

22.0 

500 


60.0 


DIAHETER (microns) 
b. LI/VPS Results 


(CASE II) 

Test Conditions: Radial Position = 2 in. 
Air Pressure = 45 psia 
Water Pressure = 115 psia 
Water Flowrate = 0.06 gal/min. 
Axial Position from Nozzle = 2 ft. 


Figure 3.37: MOD-1 Nozzle Comparison 


82 






814 


30.9 

Diane ter (nicroneters) 


Host Probable Di a= 8.1 
flrithnetic Ho an (Dig : H.l 
Area Hoan (Dig = 8.3 
Dolune Hoan (D39): 18.7 
Sauter Hoan (D32)= 24.7 


Corrector Count: 9913 


File: NAS4 36.MT 
Date: 88-94-1987 


A tup: X22S8 
lino: 99:15:26 


Total Count: 
Bun lino: 


’a: 


796 


(« 

seconds 



389 Velocity Hoan = 9.72 
VS velocity = 3.3 


si-m 


Velocity 


14.9 25.9 

(neters/socond) 


Pta= 

Ovr= 892 
llndr 235 ... 

Motr 1 Q«}: U 1 1— 1 

rital UA s 4973 BOUHSPTU 


3953 


a. P/DPA Results 



(CASE II) 

Test Conditions: Radial Position = 2 \ in. 
Air Pressure = 45 psia 
Water Pressure = 115 psia 
Water Flowrate = 0.06 gal/min. 
Axial Position from Nozzle = 2 ft. 

Figure 3.38: MOD-1 Nozzle Comparison 
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(CASE II) 

Test Conditions: Radial Position = 3 in. 
Air Pressure = 45 psia 
Water Pressure =115 psia 
Water Flowrate = 0.06 gal/min. 
Axial Position from Nozzle = 2 ft. 


Figure 3.39: MOD-1 Nozzle Comparison 


84 
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b. LI/VPS Results 


(CASE II) 

Test Conditions: Radial Position = 3^ in. 
Air Pressure = 45 psia 
Water Pressure = 115 psia 
Water Flowrate = 0.06 gal/min. 
Axial Position from Nozzle = 2 ft. 

Figure 3.40: MOD-1 Nozzle Comparison 
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To study the aforementioned results, the arithmetic mean diameter and Sauter mean diameter 
from each test were graphed as functions of radial position (Figs. 3.41 and 3.44) for each case. The 
choice of the arithmetic and Sauter mean diameters in the graphs was made to examine the count 
vs. particle size distribution. The distribution shape most associated with aerosol spray analysis 
is similar to a log-normal distribution where the distribution mode leans toward the low side of 
the distribution and conversely the distribution tail shifts to the high side of the distribution. The 
distribution is reproduced by the fact, that the arithmetic mean diameter is proportional to the 
mode of the distribution and the Sauter mean diameter is indicative of the distribution’s tail. With 
the above technique, the comparison of results from the P/DPA and the LI/VPS was performed. 

3.2.1 Discussion of Results for Comparison - CASE I 

Referring to Table 3.5, the arithmetic mean diameters measured by the LI/VPS remained approx- 
imately constant from 9.5 /im at the centerline to 10.7 fim at the edge of the spray, while the 
P/DPA values varied from 12.3 /im at the centerline to 8.8 fim at the edge of the spray. Figure 
3.41 shows the general trend in the LI/VPS and P/DPA arithmetic mean diameter to be very 
similar with a maximum deviation of 2.8 /im at the centerline and a minimum deviation of 0.1 
//m at the 2.0 inch location. Figure 3.42 shows the trend in the Sauter mean diameter to be also 
similar for both instruments. The maximum deviation is 2.2 fim at the 1.0 inch radial position 
while the minimum deviation is 0.0 for the 2.5 inch position. The maximum deviation of 2.8 fim 
for the arithmetic mean diameter, and 2.2 fim for the Sauter mean diameter can be explained as a 
result of the difference in instrument operation (automatic imaging vs light scattering and spatial 
vs temporal), the depth of field correction used by the P/DPA and no correction for the LI/VPS 
system, and to the LI/VPS instrument calibration error calculated to be ± 2.6 fim with a standard 
deviation of ± 2.0 fi m. 

3.2.2 Discussion of Results for Comparison - CASE II 

Referring to Fig. 3.43 and Table 3.5, the maximum deviation in arithmetic mean diameter of 7.1 
fim occurred at the centerline with the minimum deviation of 1.4 fim at the edge of the spray. As 
in CASE I, the LI/VPS arithmetic mean diameters remained approximately constant from 11.9 fim 
at the centerline to 12.4 fim at the edge of the spray, and the P/DPA values varied from 19.0 fim 
at the centerline to 13.8 fim at the outer edge. Figure 3.44 showed a very similar trend in Sauter 
mean diameters as a function of the radial location for both instruments. A maximum deviation 
of 6.7 fim occurred at the centerline of the spray and a minimum deviation of 0.4 fim at the 1.0 
inch location. 

In CASE II, the increase in water pressure may increase the turbulence in the outer region 
of the spray plume, which in turn caused recirculation of particles through the overlapping probe 
volumes. In addition to the explanations given in CASE I for the the differences in the arithmetic 
mean diameters we believe that since the trend for both cases is very similar (i.e., LI/VPS values 
remained approximately constant across the spray plume, while the P/DPA values decreased as 
the measurements approached the outer edge of the spray), some of the differences is due to the 
more difficult test conditions of CASE II. As we approach the outer edge of the spray, there is 
better agreement in the arithmetic mean diameter for both instruments. A possible explanation is 
the way the P/DPA operates. Recalling from Section 2.1, for proper operation of the P/DPA, the 
drops must pass through the probe volume perpendicular to Doppler fringes. Drops exactly at the 
centerline of the spray will almost always be perpendicular to these fringes and as we approach 
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the outer edge, the drops at these locations will have different directions. The result is an increase 
in run time which for CASE II varies from 2.0 sec at the centerline to 21.2 sec at the edge of 
the spray. The increase in time is an indication that more particles were rejected; therefore, the 
svstem becomes more selective and perhaps explains the smaller arithmetic mean diameter as the 
ed<*e of the spray is approached. The difference in arithmetic mean diameters in the inner region 
of the spray is attributed to the loss of small particles due to the presence of high number of liquid 
particles per volume of air which produces overlapping signals in the P/DPA. The number density 
at the center of the spray was 6970 particles/cm 3 compared to 1070 particles/cm at the edge. 
According to Dodge et al[22], by comparing the AMD with the SMD for each case, the differences 
in the shape of the distribution can be observed. Studying Figures 3.41, 3.42, 3.43, and 3.44 it is 
observed that the Sauter mean diameter compared more closely than the arithmetic mean diameter 
which suggests a difference in distribution shape for each case. 
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Test Conditions: 

Air Pressure « 65 psia. 
Water Pressure - 105 psia. 
Water Flow- rate - 0.069 gal/min. 


Figure 3.41: Comparison of Arithmetic Mean Diameters for MOD-1 Nozzle Comparison Test - 
CASE I 
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Figure 3.42; Comparison of Sauter Mean Diameters for MOD-1 Nozzle Comparison Test - CASE I 
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Figure 3.43: Comparison of Arithmetic Mean Diameters for MOD-1 Nozzle Comparison Test - 
CASE II 
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Figure 3.44: Comparison of Sauter Mean Diameters for MOD-1 Nozzle Comparison Test - CASE 
II 
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Table 3.5: MOD-1 Nozzle Comparison Results 


CASE I 
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Section 4 


CONCLUSIONS AND RECOMMENDATIONS 


This section presents the conclusions of the experimental findings and suggestions for ^ 
the experimental apparatus and drop-sizing instrumentation in future studies. The first section 
deals with the revisions to the LI/VPS, including the upgrade to dynamic particle sizing, the 
development of the calibration procedure, and the software updates. The second sect, on deals with 
the comparison of the LI/ VPS and the P/DPA, and observations concerning .the,, proper ■ 
set-up, and limitations. The final section of pertains to the improvement of the LI/VPS i to a more 
complete drop-sizing instrument, the continuation of aerosol spray analysis on the MOD-1 nozzle, 
and general observations concerning the continuing work in aerosol drop-sizing. 


4.1 LI/VPS 

The LI/ VPS has been upgraded to a system capable of performing drop-sizing analysis on dynamic 
particles. With the addition of the AD/DA converter board to the control computer, the PSP has 
shown the capability to distinguish drop-size and focus on dynamic particles m the SPM M 
freeze frame analvsis). Therefore, the LI/VPS’ drop-sizing method and focus criteria, dex eloped 
prior to this work, remains essentially intact with minor modifications. 

A two-dimensional calibration procedure for LI/VPS has been developed which allows for a 
straight-forward, step-by-step process in determining the micron/pixel correction factors associate 
with the lens magnification and camera tube non- Uneaxities With the 

procedure and the availability of a f/8 lens (he., approximate LI/VPS magnification of 1000), the 
lower-limit on the measurable, focus- dependent size-span of the LI/VPS has been reduced from 

Inchided in the LI/VPS upgrade has been the development of the PSP set-up sub-program and 
a drop-size distribution graphics display package. Due to the variability of conditions in ^rosd 
spray analysis and the flexibility of the LI/VPS, the set-up sub- program was deve loped o aid 
the operator in his decision process and allow for utilization of the full capabilities of the LI/ S. 
The add-on of the graphicW was necessary to further the Ll/VPS ab.l.ty to c ^'ze 

aerosol sprays. The graphical representation of the drop-s,ze data was ““ “ “ s 3t“e P/DPA 
specifying the proper drop-size range and a tool m the comparison of the LI/VPS and V/Ui'A. 

4.2 LI/VPS and P/DPA Comparison 

The Ll/VPS and the P/DPA compared favorably in tests performed both on the VOAG i as well 
as on the MOD-1 nozzle. Results of calibration runs performed with the \ OAG for ca.cs _»th 
and without particle dispersion showed agreement between instruments witlun ± 2.6 ,/m. 
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standard deviation of the calibration test results were all under 2.0 /xm. The small standard 
deviation indicates the accuracy of these instruments for similar test conditions. The MOD- 1 nozzle 
experiments also showed similar agreement between instruments. Results from CASE I shows a 
maximum 2.8 /xm difference in AMD and a 2.2 /xm difference in SMD. AMD values determined 
for CASE II show a higher deviation than CASE I (7.1 /jm and 2.8 /xm respectively). The AMD 
values agree quite well for the outer region of the spray where the PDPA system becomes more 
selective as explained in section 3.3.2. The SMD for both instruments follows the same general 
trend across the spray with a maximum deviation of 6.7 /xm. Considering the difference in the basic 
sizing methods employed by the two instruments and the very difficult test operating conditions, 
the LI/VPS and the P/DPA comparative measurements were surprisingly close especially for the 
SMD. 

Proper operation and set-up of the LI/ VPS and the P/DPA depend highly on the operating 
conditions specified in each test case. For this discussion, the MOD-1 nozzle is of prime interest. 
The operating conditions of the MOD-1 nozzle for the aforementioned cases, were not ideal for 
either instrument. Since the LI/VPS has limited lower size measurement capabilities, the AMD 
and SMD values determined may be slightly higher than the actual values. On the other hand, 
turbulence at the outer regions of the spray plume seemed to cause the P/DPA to reject a high 
number of counts. It is important that the operator monitor each instrument in characterizing 
any unknown aerosol spray. Even though the LI/ VPS and the P/DPA agree remarkably well, each 
instrument performs better under different test conditions. The LI/ VPS performs well in a high 
density aerosol spray, whereas the P/DPA under similar conditions, appears to have difficulties 
due to the overlap of signals (multiple particles in probe volume). Particle rejection in the P/DPA 
appears to limit the capability of this instrument to make liquid water flux measurements for the 
test conditions considered here. The P/DPA is much faster than the LI/VPS which allows for more 
versatility especially in sparse sprays. Also, the P/DPA is capable of making velocity measurements 
concurrently with the drop-size measurement, but as was shown for the MOD-1 comparison, the 
recirculation of drops associated with the turbulent spray resulted in numerous rejections. 

4.3 Suggestions and Recommendation for Future Work 

The LI/VPS, as particle sizing instrument, has progressed in stages of development. The next 
stage of development should be to upgrade the system to off-line analysis (e.g., frame storage on 
a read-write laser disk recorder), as well as increasing the program speed through hardware and 
software modifications. A study should be performed to determine the feasibility of frame storage, 
and if necessary, the error associated with such storage. The control computer, the behavior 
of imaging laser, and the PSP program structure should be studied to increase the operating 
speed of the LI/VPS. With the addition of the Micro- VAX computer, the control computer should 
not be the limiting parameter in program speed. The PSP trigger to the imaging laser doesn’t 
function consistently which makes it necessary to check for appropriate background level before 
processing. Therefore, with proper operation of the imaging trigger, unnecessary processing time 
can be avoided. Finally, to increase the speed of the LI/VPS, the PSP should be stream-lined. 
For example, the double- threshold used to determine BGL parameter for particle focus should be 
consolidated into a single threshold. 

The research on the MOD-1 nozzle and the comparison of the LI/VPS and the P/DPA should 
be continued. Operating conditions for the current work were specified by NASA. Future work on 
the MOD-1 nozzle should involve tests performed at lower water and air nozzle pressures. These 
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operating conditions would produce a larger drop-size and reduce turbulence in the spray. Also, a 
position closer to the nozzle would produce a higher number density spray which would be ideal 
for the LI/ VPS. The use of the P/DPA 200 mm transmitter lens would reduce the probe volume 
which, in turn, would reduce the probability of multiple particles in the probe volume produced 
by the high number density of droplets. The 200 mm transmitter lens was not used in the current 
work since similar sized probe volumes were needed in making the simultaneous and o\erlapping 
probe volume analysis. The above suggestions are included to improve the functionality of the two 
instruments in future studies. 

The current research and other comparison work by Dodge et al. [22] and Jackson et al. [23] 
improve the understanding of the various types of sizing techniques and assist in the development 
of accurate sizing instrumentation. The selection of a calibration/verification method or standard 
should be found for all drop-sizing instruments. The selection should be a priority for researchers 
and instrument manufacturers. 
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Section 6 


APPENDIX A: EQUIPMENT LISTING 


Device 

Manufacturer 

Model # 

Serial # 

P/DPA Transmitter 

Aerometrics Inc. 

1100 

101 

P/DPA Receiver 

Aerometrics Inc. 

2100 

101 

P/DPA Signal Processor 

Aerometrics Inc. 

PDP 3100 

103 

P/DPA Control Computer 

IBM Corp. 

AT-5170 

01619045170 

P/DPA Output Printer 

Hewlett Packard Corp. 

2225A 

261 7S3 14 1 1 

LI/ VPS Imaging Laser 

Energy Systems Inc. 

N2-50 


LI/VPS Imaging Laser 
Power Supply 

Laser Systems Inc. 

N2-50 


LI/VPS Imaging 
Laser Control Module 

Laser Holography Inc. 

N2-50 


LI/ VPS Laser 
Sync Generator 

Laser Holography Inc. 

N2-50 


LI/VPS (back-up) 
Imaging Laser 

Molectron Corp. 

UV-12 

198 

LI/ VPS (back-up) Imaging 
Laser Vacuum Pump 

Busch Inc. 

V-20 

116 

LI/ VPS Video Camera 

COHU Inc. 

2006-011 

625810 

LI/VPS Video Camera 
Control Unit 

COHU Inc. 

7910B-011 

112698 

LI/ VPS Control Computer 

Digital Equipment Corp. 

PDP11/73 


LI/VPS Output Printer 

Digital Equipment Corp. 

LA75-A2 

03555 

LI/ VPS Video Processor 

Recognition Concepts Inc. 

TRAPIX 55/32Q 

134 

Computer Terminal 

Digital Equipment Corp. 

VT-240 

HK14705 

Computer Terminal 

CIE Terminals Inc. 

CIT-220+ 

8471C6916 
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Video Monitor 

Video Monitor 

Video Cassette Recorder 

Video Cassette Recorder 

Laser/optical 

Disk Recorder 

Decwriter 

Digital Oscilloscope 
Measurement 
Plotting System 
Real-time Oscilloscope 
Digital Multimeter 
Pressure Transducer 
Pressure Transducer 
Direct Current 
Power Supply 
Flowmeter-regulator 
Monodisperse 
Drop Generator 
Test Nozzle 
Air and Water Supply 
Pressure Regulator 
Water Pressure Vessel 
Isolation Table 
High-pressure 
Charging Hose 
Dead Weight 
Pressure Tester 
Balance Scale 


Sanyo Corp. 

SONY Corp. 

RCA Corp. 

Panasonic Corp. 

Panasonic Corp. 

Digital Equipment Corp. 
Hewlett-Packard Corp. 
Hewlett-Packard Corp. 

Tektronix Inc. 

John Fluke Mfg. Co., Inc. 
OMEGA Eng., Inc. 
OMEGA Eng., Inc. 
Hewlett-Packard Corp. 

Cole-Palmer Inc. 

TSI Inc. 

NASA Lewis Space Center 
University of Nebraska 
WATTS Regulator Co. 
Brunswick Corp. 

Newport Research Corp. 
Yellow Jacket Inc. 

American Steam Gage Co. 

HOWE Inc. 


AVM255 

55805757 

CKV-1900F 

204071 

VET650 

1032FM243 

NV-8950 

B5HL00491 

TQ-2320F(A) 

EH4669001 

LA120AA 

PNE1366 

54200A 

2511A-00639 

7090A 

2430A00344 

549 

7365 

8024 B 

3715516 

PX304-150 

850502 

PX304-150 

850311 

6200 B 

24 11 A- 12365 

PR004 

FM044-40G 

345000 

167 

MOD-1 


2-26A 

8305 

6121020-024 

1003 

WPP0031A 



3074131 
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Section 7 


APPENDIX B Design and Implementation of the PSP Laser Trigger 


Due to the availability of the existing laser sync circuit (LSC) and the AD/DA converter board, 
the development of the PSP software generated trigger was simplified. With the aforementioned 
hardware, the PSP software, utilizing available FORTRAN callable commands, directs a digital 
value to the AD/DA board. The AD/DA board converts the digital value to the appropriate analog 
signal. The analog signal is then sent to the LSC. The analog signal from the control computer is 
paralleled with the sync signal from the CCU, and the resulting signal triggers the imaging laser. 
The above process was used as the basis for LI/VPS conversion from the CPM to the SPM. Except 
for cabling, the majority of work in the modification dealt with the LSC. Figure B7.1 shows the 
overall circuitry of the LSC with special attention given to the source of the LSC laser trigger and 
the position marked by the Xs. The major addition to the LSC circuitry was the two AND gates 
(Appendix B, Fig. B7.2. Therefore an analog signal from the control computer must be present 
at the first AND gate before the imaging laser can be triggered. The above method, therefore, 
facilitates SPM operation for the LI/VPS. 
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Section 8 


APPENDIX C.l: PSP Set-up Program 


C 


PROGRAM MENU 


+ 4 ++++++++++ C 

C PROGRAM DEVELOPED TO SET-UP OPERATING PARAMETERS FOR THE 
C THE PARTICLE SIZING PROGRAM UTILIZING A MENU-TYPE FORMAT. 
C 


44444444444 C 

REAL*8 PROCESS (2) , YESN0(2) , ADVANCE(2) , LIMIT(3), 
B0UNDRY(2) REAL*8 MAG(2) , A, B, C, D, F, G, H, I, J, V, W 

INTEGER TOD (4), D0Y(5) 

CHARACTER* 10 NUMBER 
CHARACTER+6 FILE1 , FILE2 


CHARACTER* 1 L(6) , M(6) 

LOGICAL* 1 IKEY 
BYTE ESC, LINE(50,4) 

DATA ESC / 27 / 

DATA PROCESS, YESNO / ’STATIC’, ’ DYNAMC * , ’YES ’, ’NO 


’ / DATA ADVANCE, FILE1 / 

DATA LIMIT, FILE2 / ’TIME ’ 
DATA BOUNDRY, MAG / ’PROCSS 
DATA IGPST, WIDTH, NGRPS, 
2000 / DATA JXSTR, JXDST, 

DATA A, B, C, D / ’DYNAMC 
DATA F, G, H, I / ’REJECT 
DATA J, V, W / ’NO ’, 
DATA LINE / 200*’ ’ / 

DATA NUMBER / ’1234567890 
EQUIVALENCES FILE1) 
EQUIVALENCES FILE2) 

CALL DATE(DOY) 

CALL TIME(TOD) 

WRITES 1) ESC, ESC, ESC, 


’SINGLE’, ’AUTO ’, ’TEMP01’ / 

’FRAME ’, ’ PARTCL ’ , ’ TEMPO 1’ / 

, ’REJECT’, ’LOW ’, ’HIGH ’ / 

ITHRSH, LIMVAL / 5, 5.0, 50, 90, 
JYSTR , JYDST / 50, 400, 50, 400 / 
, ’YES ’, ’AUTO ’, ’PARTCL’ / 

, ’YES ’, ’YES ’, ’YES ’ / 

LOW ’ , ’NO ’ / 

/ 


ESC 


103 



1 F0RMAT(1X, A1 , ’ [81’ ,A1, ’ [?251’ ,A1 , ’ [H’ ,A1, ’ [2J’> 
C 


C*** IF DATE NOT SET, INSERT "NO DATE" INTO DATE FIELD 
C 


C 

C 

C. 

C 


IF ( DOY(l) .EQ. 
IF( DOY(l) .EQ. 
IF( DOY(l) .EQ. 
IF ( DOY(l) .EQ. 
IF ( DOY(l) .EQ. 


’00’ ) DOY (5) = ’ ’ 
’00’ ) DOY (4) = ’TE’ 
’00’ ) DOY (3) = ’DA’ 
’00’ ) DOY (2) = ’0 ’ 
’00’ ) DOY(l) = ’ N’ 


INPUT PREVIOUSLY STORED PSP PARAMETERS. 


0PEN(UNIT=1, FILE=’SETUP.MNU’ .STATUS* ’OLD’ ,ERR= 5) 
READ(l.lO) A. B, C, D, LIMVAL, F, G, H, FILE1 , I, J, FILE2, 
t IGPST , WIDTH, NGRPS , JXSTR, JXDST, JYSTR , JYDST , ITHRSH, 

V, W READ(1 ,2) 

2 F0RMAT(2(/)) 

DO 4 11=1,4 

READ(1 ,3) (LINE(JJ.II), JJ=1,50) 

3 FORMAT (3X.50A1) 

4 CONTINUE 

GOTO 9 5 OPEN (UNIT=1 , FILE= ’SETUP. MNU’ , STATUS* ’ NEW ’ ) 
WRITE(l.lO) A, B, C, D, LIMVAL, F, G, H, FILE1 , I, J, FILE2, 
ft IGPST, WIDTH, NGRPS, JXSTR, JXDST, JYSTR, JYDST, ITHRSH, V, W 
WRITE (1,2) 

DO 6 11=1,4 

WRITECl ,3) (LINE(JJ.II), JJ=1,50) 

6 CONTINUE 

9 WRITE(7,10) A, B, C, D, LIMVAL, F, G, H, FILE1 , I, J, 

FILE2, ft IGPST, WIDTH, NGRPS, JXSTR, JXDST, JYSTR, JYDST, 

ITHRSH, V, W CL0SE(UNIT=1) 

10 FORMATC ’+’,T39,’ SETUP ’/T26,’ PARTICLE SIZING PROGRAM (ver. 

4)’// ft 33( ’ - ’ ) , ’ PROCESSING OPTIONS’ ,T54,2S(’-’)/ 

* T3, ’ (A) ’ ,T9,A6,T19, ’Type of Processing* ,T49, 
ft ’(STATIC/DYNAMIC)’,/ 
ft T3, ’ (B) ’, T9, A6.T19, ’Focus Criteria* ,T49, 
ft ’(YES/NO)’,/ 

ft T3, ’ (C) ’ ,T9,A6,T19, ’Type of Frame Advance* ,T49, 
ft ’(AUTO/SINGLE)',/ 

ft T3, ’ (D) ’, T9, A6.T19, ’Processing Limit *,T49, 
ft ’ (TIME/FRAME/PARTICLE) ’ , / 

ft T3,’(E) ( ’ ,T9,I6,T15, ’) Limiting Value’ ,T49, 
ft ’ (seconds/f rames/particles) ’ , / 
ft T3, ’ (F) ’ ,T9,A6,T19, ’Boundary Particles* ,T49, 
ft ’ (PROCESS/REJECT) ’ ,/ 
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t 36( ’-’),’ OUTPUT OPTIONS ’ ,26(’-’)/ 

* T 3 , ’ (G) ’, T9, A6.T19, ’General Results (to PRINTER) 
(YES/NO)’/ * T19, ’WRITE TO FILE (YES/NO) ’ ,T49, 

ft »(K) FILE HEADER (4 lines)’,/ 

ft T3, ’ (H) ’, T9, A6.T19, 'Average Particle size data -- ’,T49, 
ft ’(L) FILE: (’ ,A6,’).0UT’ ,/ 

ft T3, ’ (I) ’ ,T9,A6,T19, ’Group Breakdown data ’ ,9( ’-’),’/’ / 
ft T3 , ’ ( J) ’ ,T9 , A6.T19 , ’ Per Frame data ’, 13( ’-’),’>’ ,T49 , ft 

’(M) FILE: (’,A6,’) .DAT’,/ 

ft 35( ’ - ’ ) , ’ GENERAL OPTIONS ’,26(’-’)/ 
ft T3,’(N) Group Start =( ’,13,’) ’ 
ft T29 , ’ (0) Group Width =(’,F4.1,’) ’ 
ft T55, ’ (P) # of Groups = ( ’ ,13, *) */ 
ft T3,’(Q) X Window Start * (’,13,*)’, 
ft T49 , ’ (R) X Window Width - (’,13,’)’/ 
ft T3,’(S) Y Window Start * (’,13,’)’, 
ft T49 , ’ (T) Y Window Width = (’,13,’)’/ 
ft T3,’(U) Threshold = (’,13,’)’ 
ft T29, ’ (V) Lens = \A6, 

ft T49, ’ (W) Markers = ’ ,A6, ’ (YES/NO) ’/78(’-’)) 

WRITE (7, 20) ESC, ESC 

20 FORMAT ( ’ + ’ ,A1 , ’ [23 ; 2H ’ , A1 , ’ [OJ ’ , ’ (X) to exit SETUP menu or 

ft (Z) to begin Particle Sizing Program */ 

ft T3, ’Enter Letter to change specific Parameter ??’) 


C 

C 

c 

c 

30 


SPECIFY PSP PARAMETER W/ KEY TOGGLE OR KEYBOARD ENTRY . 


CALL IP0KE("44," 10000 . 
IKEY = ITTINRO 

OR. 

IF(IKEY.LT.O) 

GOTO 

30 

IF(IKEY.EQ.’A’) 

GOTO 

100 

IF(IKEY.EQ.’B’) 

GOTO 

200 

IF(IKEY.Eq.’C’) 

GOTO 

300 

IF(IKEY.EQ.’D’) 

GOTO 

400 

IF(IKEY.EQ.’E’) 

GOTO 

500 

IF(IKEY.EQ.’F’) 

GOTO 

600 

IF(IKEY.EQ.’G’) 

GOTO 

700 

IF(IKEY.EQ.’H’) 

GOTO 

800 

IF (IKEY. EQ. ’I’) 

GOTO 

900 

IF(IKEY.EQ. ’ J’) 

GOTO 

1000 

IF(IKEY.EQ.’K’) 

GOTO 

1100 

IF(IKEY.EQ.’L’) 

GOTO 

1200 

IF(IKEY.EQ.’M’) 

GOTO 

1300 

IF(IKEY.EQ.’N’) 

GOTO 

1400 

IF(IKEY.EQ.’O’) 

GOTO 

1500 


IPEEK("44) ) 
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IF(IKEY.Eq.’P’) GOTO 1600 

IFUKEY.EQ.’CT) GOTO 1700 

IF(IKEY.EQ.'R’) GOTO 1800 

IF(IKEY.EQ.’S’) GOTO 1900 

IF(IKEY.EQ.'T’) GOTO 2000 

IF(IKEY.Eq.'U’) GOTO 2100 

IF(IKEY.Eq.’V’) GOTO 2200 

IF(IKEY.Eq.’W’) GOTO 2400 

IF(IKEY.Eq.’X’) GOTO 2600 

C IF(IKEY.Eq. ’Y*) GOTO 2500 

IF(IKEY.Eq.’Z’) GOTO 2600 

GO TO 30 

90 FORMAT ( ’ + ’ ,11, » [23;2H* ,11. ‘ [OJ’ , 

t ’Enter new value for (*,11,*) here **>’,$) 
C 

C... (A) SPECIFY PROCESS TYPE (DYHAMIC/STATIC) 

C 

100 CALL IP0KE("44, "167777 .AND. IPEEK("44)) 
IF(A.Eq. PROCESS (1)) GOTO 150 
A = PROCESS d) 

GOTO 180 

150 A = PROCESS (2) 

180 WRITE(7 , 190) ESC, A 
190 FORMAT ( ’ + ’ , A 1 , ’ [5 ; 8H ’ , A6) 

GOTO 30 C 

C... (B) SPECIFY FOCUS (YES/NO) 

C 

200 CALL IPOKE( "44, "167777 .AND. IPEEK("44)) 
IF(B.Eq.YESNOd)) GOTO 250 
B - YESNO(l) 

GOTO 280 

250 B * YESN0(2) 

280 WRITE (7, 290) ESC, B 

290 F0RMAT(’+’,A1,’[6;8H’,A6) 

GOTO 30 

C 

C... (C) TYPE OF FRAME ADVANCE (SINGLE/ AUTO) 

C 

300 CALL IP0KE("44," 167777 .AND. IPEEK("44)) 

IF (C . EQ . ADVANCE( 1) ) GOTO 350 
C = ADVANCE(l) 

GOTO 380 

350 C = ADVANCE (2) 

380 WRITE (7, 390) ESC, C 

390 FORMAT ( ’ + ’ ,A1 , ’ [7 ;8H ’ , A6) 

GOTO 30 


106 


c... (D) PROCESSING LIMIT (TIME/PARTICLE/FRAME) 

C... NOTE: DUE TO COMPUTER LIMITATIONS TIME IS NOT INCLUDED C 
400 CALL IP0KE("44," 167777 .AND. IPEEKO'44)) 

IF (D . EQ . LIMITC 1) . OR . D . EQ . LIMIT (2) ) GOTO 450 
D = LIMIT(l) 

GOTO 480 

450 IF(D . EQ .LIMIT (2) ) GOTO 470 
D = LIMIT(2) 

GOTO 480 

470 D = LIMITC3) 

480 WRITE (7, 490) ESC, D 
490 FORMATC ,A1 , ’ [8;8H’ ,A6) 

GOTO 30 

C 

c... (E) SPECIFY LIMITING VALUE 

C 

500 CALL IPOKE ("44, "167777 .AND. IPEEK("44)) 

WRITE(7 ,90) ESC, ESC, ’E’ 

READ (5,*) LIMVAL 
580 WRITE (7, 590) ESC, LIMVAL 
590 FORMAT( ’+’ ,A1 , ’ [9;8H’ ,16) 

WRITE(7 ,20) ESC, ESC 
GOTO 30 
C 

C. . . (F) BOUNDARY ANALYSIS 

C... NOTE: UNAVAILABLE 
C 

600 CALL IPOKEC "44, "167777 .AND. IPEEK("44)) 
IF(F.EQ.BOUNDRYCl)) GOTO 650 
F <= BOUNDRY(l) 

GOTO 680 

650 F = B0UNDRYC2) 

680 WRITEC7.690) ESC, F 

690 FORMATC ,A1, ’ [10;8H’ ,A6) 

GOTO 30 
C 

C. . . CG) OUTPUT GENERAL RESULTS TO PRINTER 
C 

700 CALL IP0KEC"44," 167777 .AND. IPEEKC"44)) 
IFCG.EQ.YESNOCl)) GOTO 750 
G = YESNOCl) 

GOTO 780 

750 G = YESN0C2) 

780 WRITEC7 ,790) ESC, G 

790 FORMATC ,A1 , ’ [12; 8H’ , A6) 
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GOTO 30 


C 

C... (H) WRITE TO FILE: ANALYSIS SUMMARY (YES/NO) 

C 

800 CALL IP0KEC44, "167777 .AND. IPEEK("44)) 

IF(H .EQ . YESNO(l) ) GOTO 850 
H = YESNO(l) 

GOTO 880 

850 H = YESN0(2) 

880 WRITE (7, 890) ESC, H 

890 FORMAT ( ’ +’ ,A1 , ’ [14;8H’ ,A6) 

GOTO 30 
C 

C... (I) WRITE TO FILE: GROUP BREAKDOWN (YES/NO) 

C 

900 CALL IPOKEC "44, "167777 .AND. IPEEK("44)) 

IF (I .EQ . YESNO(l) ) GOTO 950 
I = YESNO(l) 

GOTO 980 

950 I « YESN0C2) 

980 WRITE (7, 990) ESC. I 

990 F0RMATC’+’,A1,’[15;8H’,A6) 

GOTO 30 
C 

C... (J) WRITE TO FILE: PER FRAME DATA (YES/NO) 

C 

1000 CALL IPOKEC "44, "167777 .AND. IPEEK("44)) 

I F ( J . EQ . YESN 0(1)) GOTO 1050 
J = YESNO(l) 

GOTO 1080 
1050 J = YESN0C2) 

1080 WRITE(7 , 1090) ESC, J 
1090 FORMAT ( ’ + ’ ,A1 , ’ [16;8H’ ,A6) 

GOTO 30 
C 

C... (K) SPECIFY FILE HEADER 

C 

1100 CALL IPOKEC "44, "167777 .AND. IPEEKC"44)) 

WRITEC7, 1110) ESC, ESC, ESC, ESC 

1110 FORMAT C IX, A1 , ’ [81* ,A1 , ’ [7251* ,A1 , ’ [H* ,A1, ’ [2J* ) 
WRITEC7 , 1111) CCLINECJJ.II), JJ=1,50), 11=1,4) 

1111 F0RMATO + ’ ,T39, ’SETUP ’/T35, ’FOR FILE HEADER’// 

ft T36, ’CHANGE CY/N) ’ ///T16 ,54C ’ * ’ ) ,/T16, ,T69, ’*’ , 
ft 4C/.T16,’* ’ 50A1 , ’ * ’ ) ,/T16 , ’ * ’ ,T69, , 
ft /T16,54C’*’)//T31,’C4 LINES/50 SPACES each)’) 

CALL IP0KEC"44," 10000 .OR. IPEEKC"44)) 
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1115 IKEY » ITTINRO 

IFCIKEY.LT. 0) GOTO 1115 
IFCIKEY.EQ.’Y’) GOTO 1120 
IF(IKEY.EQ.’N’) GOTO 1190 
GOTO 1115 

1120 CALL IP0KE("44," 167777 .AND. IPEEK("44)) 

DO 1140 11=1,4 

WRITE (7, 1125) ESC, ESC, II 

1125 FORMAT ( ’ +’,A1,’[17;3H’,A1,’ [OJ * , 
jt ’Change Line C’,11,’), (Y/N)’) 

CALL IP0KE("44," 10000 .OR. IPEEK("44)) 

1126 IKEY = ITTINRO 
IFCIKEY.LT. 0) GOTO 1126 
IFCIKEY.EQ.’Y’) GOTO 1127 
IFCIKEY . EQ . ’ N ’ ) GOTO 1140 
GOTO 1126 

1127 CALL IP0KEC n 44, "167777 .AND. IPEEKC"44)) 
WRITEC7 , 1130) ESC, ESC, ESC, II 

1130 FORMATC * + ’ , A1 , ’ [22 ; 1H ’ , A1 , ’ [0 J ’ , A1 , ’ [?8h’ , 
* ’ Line (Ml,’) —>’,$) 

IF(II.EQ.l) GOTO 1141 
IFCII.EQ.2) GOTO 1143 
IFCII.EQ.3) GOTO 1145 
IFCII.EQ.4) GOTO 1147 

1140 CONTINUE 
GOTO 1180 

1141 READC5.1170) (LINE(JJ.l), JJ=1,50) 

WRITEC7, 1142) ESC, (LINE(JJ.l), JJ=1.50) 

1142 FORMATC ’+’ ,A1,’ [9; 15H’ ,’* ’,50A1,’ *’) 
GOTO 1140 

1143 READ(5 ,1170) (LINECJJ.2). JJ-1,50) 

WRITE(7 , 1144) ESC, (LINE(JJ,2), JJ=1,50) 

1144 FORMATC ’+’,A1,’[10; 15H’ , ’* ’,50A1,’ *’) 
GOTO 1140 

1145 READ(5 ,1170) (LINE(JJ,3), JJ=1,50) 

WRITEC7 , 1146) ESC, (LINE(JJ,3), JJ=1,50) 

1146 FORMATC ’+’,A1,’[11;15H’,’* *’) 

GOTO 1140 

1147 READ(5 ,1170) (LINE(JJ,4), JJ-1,50) 

WRITE(7 ,1148) ESC, (LINE(JJ,4), JJ=1,50) 

1148 FORMATC ’+’ ,A1, ’ [12; 15H’ ’.SOAl,’ *’) 
GOTO 1140 

1170 FORMAT (50A1) 

1180 GOTO 1100 

1190 WRITEC7.1) ESC, ESC, ESC, ESC 
GOTO 9 
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c 

C... (L) FILE SPECIFICATION: GENERAL ft GROUP DATA 

C 

1200 CALL IPOKE( "44, "167777 .AND. IPEEK("44)) 
WRITE(7, 1220) ESC, ESC, 'L' 

1220 FORMATC ’+’ ,A1 , ’ [23; 2H’ , A1 , ’ [OJ ’ , 

* ’File Name C’.Al,’) (4 letters) here ==>’,$) 
READ (5, 1230) (L(II), 11=1,4) 

1230 F0RMAT(4A1) 

WRITE(7 , 1240) ESC, (L(II), 11=1,4) 

1240 FORMAT(’+’ ,A1, ’ [14;60H* ,4A1) 

WRITEC7.1250) ESC, ESC, 'L' 

1250 FORMATC ♦ ’ ,A1, ’ [23;2H» ,A1,» [OJ' , 

ft 'File Number (’,A1,') (2 numbers) here ==>’,$) 
READ(5, 1260) LI 
1260 F0RMAT(I2) 

II = Ll/10 
10 = L1-I1*10 
IFCL1.LT. 10) 11=10 

IFCIO .EQ. 0) 10=10 
L(5) « NUMBER(I1 : II ) 

L(6) = NUMBER(IO:IO) 

IFCL1.LT. 10) GOTO 1280 
WRITEC7, 1270) ESC, LI 
1270 FORMATC'+’ ,A1, ’ [14;64H* ,12) 

WRITEC7.20) ESC, ESC 
GOTO 30 

1280 WRITEC7 , 1290) ESC, LI 
1290 FORMATC ’+’ ,A1, ’ [14;64H0* ,11) 

WRITEC7.20) ESC, ESC 
GOTO 30 
C 

C... CM) FILE SPECIFICATION: PER FRAME DATA 

C 

1300 CALL IP0KEC"44, "167777 .AND. IPEEKC"44)) 

WRITEC7 , 1320) ESC, ESC, ’M’ 

1320 FORMATC ’+’ ,A1,» [23;2H* , A1 , * [OJ ’ , 

ft 'File Name C'.Al,’) C4 letters) here ==>',$) 
READ C5, 1330) (M(II), 11=1,4) 

1330 F0RMATC4A1) 

WRITEC7.1340) ESC, CMCII), 11=1,4) 

1340 FORMAT C ’ + ’ ,A1, ’ [16;60H' ,4A1) 

WRITE C7, 1350) ESC, ESC, 'M' 

1350 FORMATC ’+’ ,A1, ’ [23;2H* ,A1,’ [OJ’ , 

ft 'File Number C’.Al,’) C2 numbers) here ==>’,$) 
READ C5, 1360) Ml 
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1360 


1370 


1380 

1390 


C 

C. . . 
C 

1400 


1480 

1490 


C 

C. . . 
C 

1500 


1580 

1590 


C 

C. . . 
C 

1600 


1680 

1690 


F0RMAT(I2) 

II = Ml/10 
10 = M1-I1*10 
IF(Ml.LT.lO) 11=10 
IFCIO .EQ. 0) 10=10 
M(5) = NUMBER(I1 : II) 

M(6) = MUMBER(IO : 10) 

IF(Ml.LT.lO) GOTO 1380 
WRITE(7 , 1370) ESC, Ml 
FORMAT( ’+ ’ ,A1 , ’ [16 ; 64H ’ ,12) 

WRITE(7 ,20) ESC, ESC 
GOTO 30 

WRITE(7 , 1390) ESC, Ml 
FORMAT( ’+’ ,A1, ’ [16;64H0’ ,11) 

WRITE(7 ,20) ESC, ESC 
GOTO 30 

(N) DROP-SIZE GROUP BREAKDOWN: STARTING VALUE 

CALL IP0KE("44," 167777 .AND. IPEEK("44)) 
WRITE(7 ,90) ESC, ESC, ’N’ 

READ (5,*) IGPST 

WRITE (7, 1490) ESC, IGPST 

FORMAT( ’ + ’ ,A1 , * [18 ;21H’ ,13) 

WRITE(7 ,20) ESC, ESC 
GOTO 30 

(O) DROP-SIZE GROUP BREAKDOWN: INTERVAL WIDTH 

CALL IP0KE("44," 167777 .AND. IPEEK("44)) 

WRITE (7, 90) ESC, ESC, ’0’ 

READ (5,*) WIDTH 

WRITE(7 , 1590) ESC, WIDTH 

FORMAT ( ’ + ’ ,A1,* [18;46H* ,F4.1) 

WRITE (7, 20) ESC, ESC 
GOTO 30 

(P) DROP-SIZE GROUP BREAKDOWN: # OF GROUPS 

CALL IP0KE("44, "167777 .AND. IPEEK("44)) 
WRITE (7, 90) ESC, ESC, >P> 

READ (5,*) NGRPS 

WRITE(7 , 1690) ESC, NGRPS 

FORMAT(’+’ ,A1 , ’ [18;73H’ ,13) 

WRITE(7 ,20) ESC, ESC 
GOTO 30 
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(Q) SIZING WINDOW (PIXEL SPEC): X STARTING VALUE 


C 

C. . . 

C 

1700 CALL IPOKE( "44, "167777 .AND. IPEEK("44)) 

WRITE(7 ,90) ESC, ESC, ’Q’ 

READ (5,*) JXSTR 
1780 WRITEC7 , 1790) ESC, JXSTR 
1790 FORMAT(’+* ,A1, ’ [19;24H’ ,13) 

WRITEC7.20) ESC, ESC 
GOTO 30 
C 

C... (R) SIZING WINDOW (PIXEL SPEC): X SCREEN WIDTH 

C 

1800 CALL IP0KE("44," 167777 .AND. IPEEK("44)) 
WRITE(7,90) ESC, ESC, »R> 

READ (5,*) JXDST 
1880 WRITE(7 , 1890) ESC, JXDST 
1890 FORMAT( ’ +’ ,A1 , ' [19 ;72H’ , 13) 

WRITE (7, 20) ESC, ESC 
GOTO 30 
C 

C... (S) SIZING WINDOW (PIXEL SPEC): Y STARTING VALUE 

C 

1900 CALL IP0KE("44, "167777 .AND. IPEEK("44)) 
WRITE(7,90) ESC, ESC, ’S’ 

READ (5,*) JYSTR 
1980 WRITE(7 , 1990) ESC, JYSTR 
1990 FORMAT ( ’+ ’ ,A1 , ’ [20;26H’ , 13) 

WRITE (7, 20) ESC, ESC 
GOTO 30 
C 

C... (T) SIZING WINDOW (PIXEL SPEC): Y SCREEN WIDTH 

C 

2000 CALL IPOKE( "44, "167777 .AND. IPEEK("44)) 

WRITE (7, 90) ESC, ESC, »T» 

READ (5,*) JYDST 
2080 WRITE(7 ,2090) ESC, JYDST 
2090 FORMAT ( ’ + ’ ,A1 , ’ [20;70H* , 13) 

WRITE (7, 20) ESC, ESC 
GOTO 30 
C 

C... (U) INPUT SIZING THRESHOLD 

C 

2100 CALL IPOKE( "44, "167777 .AND. IPEEK("44)) 

WRITE(7 ,90) ESC, ESC, 'U* 

READ(5 , *) ITHRSH 
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2180 WRITE(7 ,2190) ESC, ITHRSH 
2190 F0RMAT(’ + \A1,’[21;19HM3) 

WRITE(7,20) ESC, ESC 
GOTO 30 
C 

C... (V) INPUT SYSTEM MAGNIFICATION (HIGH/LOW) 

C 

2200 CALL IP0KE("44," 167777 .AND. IPEEK("44)) 

IF(V.EQ .MAG(l)) GOTO 2250 
V = MAG(l) 

GOTO 2280 
2250 V = MAG(2) 

2280 WRITE(7 ,2290) ESC, V 
2290 FORMAT( ’ + ’ ,A1 , ’ [21 ;39H’ ,A6) 

GOTO 30 
C 

c... (W) DIAGNOSTIC MARKERS PLACED ON COUNTED PARTICLES 

C 

2300 CALL IP0KE("44, "167777 .AND. IPEEK("44)) 

IF (W . EO . YESNO ( 1 ) ) GOTO 2350 
W = YESNO (1) 

GOTO 2380 

2350 W = YESNO (2) 

2380 WRITE(7 ,2390) ESC, W 
2390 FORMAT ( ’+’ ,A1 , ’ [21 ;62H’ ,A6) 

GOTO 30 
C 

C. . . (Z) STORE SET-UP PARAMETERS AND START PSP 

C 

2600 CALL IP0KE("44," 167777 .AND. IPEEK("44)) 

0PEN(UNIT=1, FILE=’ SETUP. MNU’ , STATUS*’ NEW’) 
WRITE(l.lO) A, B, C, D, LIMVAL, F, G, H, FILE1 , I, 
FILE2, * IGPST, WIDTH, NGRPS , JXSTR, JXDST, JYSTR , 

ITHRSH, V, W WRITE(1 ,2605) 

2605 F0RMAT(2(/)) 

DO 2650 11=1,4 

WRITE(1 ,2610) (LINE(JJ,II), JJ=1,50) 

2610 FORMAT (3X.50A1) 

2650 CONTINUE 

CL0SE(UNIT=1) 

WRITE (7, 2690) ESC, ESC, ESC, ESC 
2690 FORMAT ( IX, A1 , ’ [8h’ ,A1, ’ [?25h’ ,A1, ’ [H* ,A1, ’ [2J’) 
IF(IKEY.EQ.'Z') CALL SETCMD ( * RUN PSP1’) 

CALL EXIT 
2700 STOP 
END 


J, 

JYDST, 
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Section 9 


APPENDIX C.2: PSP Graphical Presentation of Results 


C 


PROGRAM GRAPH 



++++♦+♦++♦+ C 


C PROGRAM DEVELOPED FOR THE PSP TO GRAPHICALLY REPRESENT C 

THE GROUP BREAK-DOWN DATA ON A DEC COMPATIBLE TERMINAL. C 


++++++♦♦+++ c 

DIMENSION X(IOOO) , Y(IOOO), Yl(70), Xl(70) 

BYTE ESC, CSI, TIM(9), DAY(9) 

CHARACTER* 1 A, B, C, H, II, NUMXO, NUMYO, NUMX1, NUMY1 
CHARACTER* 11 CO 

CHARACTER* 1 NAME (7) , JUNK 
CHARACTER* 10 NUMBER, FILE 
DATA A, ESC, CSI /’*’, 27, 155 / 

DATA B, C, H /’[’, ’f’/ 

DATA NUMBER / ’1234567890' / 

DATA CO, PI / ’ », 3.1415926 / 

EQUIVALENCE (NAME(l) ,B) 

EQUIVALENCE (NAME(2) ,NUMX1) 

EQUIVALENCE (NAME (3) .NUMXO) 

EQUIVALENCE (NAME(4),C) 

EQUIVALENCE (NAME(5) ,NUMY1) 

EQUIVALENCE (NAME(6) , NUMYO) 

EQUIVALENCE (NAME(7),H) 

1 DO 5 1=1,1000 

X(I)=0.0 
Y(I)=0.0 
5 CONTINUE 

DAVSUM =0.0 
DSSUM =0.0 
DVSUM = 0.0 
DWSUM = 0.0 
SUMN = 0.0 
SUM = 0.0 
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WEIGHT =0.0 
COUNT =0.0 
DO 6 1=1,70 
X1(I)=0.0 
Y1(I)=0.0 
6 CONTINUE 

C 

C... INPUT DATA FILE NAME. 

C 

WRITE(7, 10) 

10 FORMAT (IX. ’INPUT DATA FILE ==> ’,$) 

READ (5,20, ERR-200 ) FILE 
20 FORMAT(AIO) 

OPEN (UNIT=2 , NAME=FILE , STATUS- ’ OLD ' ) 

DO 30 1=1,1000 

READ (2.*, END-40) X(I),Y(I) 

30 CONTINUE 

40 INDEX = 1-1 

IFCINDEX.EQ.68) GO TO 94 

YMAX = -1.0E+30 

XMIN = 1.0E+30 

XMAX = -1.0E+30 

DO 50 K=l, INDEX 

XMIN = AMIN1(XMIN,Y(K)) 

50 XMAX = AMAX1 (XMAX.Y(K)) 

XCEN = (XMIN+XMAX)/2.0 
XSCALE = (XMAX-XMIN)/67.0 
DO 60 J=1 ,67 

X 1 ( J ) -XMIN+XSC ALE* FLOAT ( J) 

60 CONTINUE 

DO 90 L=l, INDEX 
DO 70 M=1 ,67 

IF(Y(L) .GT.Xl(M)) GO TO 70 
GO TO 80 
70 CONTINUE 

M=68 

80 Y1(M)=Y1(M)+1.0 

90 CONTINUE 

C 

C... STATISTICAL MEAN DIAMETERS DETERMINED 

C 

DO 93 1=1,67 

IF(I.EQ.l) GRPAVG = Xl(I)/2.0 
IF(I.NE.l) GRPAVG = (Xl(I)+Xl(I-l))/2.0 
DIAMAX = YMAX 

DAVSUM = DAVSUM + Y1(I)*GRPAVG 
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DSSUM 

= DSSUM 


Y1(I)*GRPAVG**2 

DVSUM 

- DVSUM 


Y1(I)*GRPAVG**3 

DWSUM 

« DWSUM 


Y1(I)*GRPAVG**4 

SUMN 

» SUMN 

+ 

Y1(I) 

SUM 

= SUM 

+ 

Y1 (I) 


WEIGHT « WEIGHT + (4 . /3 . *PI*(GRPAVG/2 . )**3)*Y1(I) 

93 CONTINUE 

DAV = DAVSUM/SUMN 
DS = SQRT(DSSUM/SUM) 

DV = (DVSUM/SUMN)** ( 1 . /3 . ) 

DVS - DVSUM/DSSUM 
DW = DWSUM/DVSUM 
NSUM * IFIX(SUM) 

GOTO 97 

94 REWIND (2) 

READ (2,*) GCD , DAV , DS , DV ,DVS , DW 
DO 95 K*1 ,67 
X1(K) - X(K+1) 

Y1(K) » Y(K+1) 

95 CONTINUE 
XMAX » Xl(67) 

XMIN = XI (1)- (XI (67) -XI (66)) 

XCEN = (XHIN+XMAX) /2 . 0 
XQUA = (XCEN-XHIN) /2 . 0 
X14 = XMIN + XQUA 
X34 = XCEN + XQUA 
C0UNT=0 . 0 

97 DO 100 K=1 ,67 

COUNT=COUNT + Y1(K) 

YMAX = AMAX 1 ( YMAX , Y1 (K) ) 

100 CONTINUE 

NSUM = IFIX(COUNT) 

YSCALE = YMAX/21 .0 

YSCAL1 * YMAX/199.0 

IYMAX = IFIX(YMAX) 

IYCEN = IFIX(.50*YMAX) 

IY14 « IYCEN/2 
IY34 * IYCEN + IY14 
C 

C... UTILIZING THE PSEUDO-GRAPHIC CAPABILITIES OF 
C... DEC COMPATIBLE TERMINALS: INITIATE GRID. 

C 

CALL CHARGR 
WRITE(7 , 105) ESC, ESC 
105 FORMATC + ’ ,A1, ’*0’ ,A1, 'n’) 

WRITE(7 , 110) ESC, ESC, ESC, ESC 
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110 


FORMAT (IX, A 1 . ’ [H’ ,A1, ’ [2J’ ,A1 , ’ [?3h» ,A1 . * C?251») 
WRITE (7, 115) C0(1 : 1) , IYMAX 

FORMATC *♦’ , 1X,A1 ,2X ,I4,2X, ' w* , 16( *<1* ) * » 16( ’<!* ) » * 

t 16(’q’) , ‘v\16(’q’).’k’) 

117 F0RMAT(2X,A1 , 2X.I4.2J, ’n’ , 16( ’q’ ) , ’n’ ,16( ’q’ ) , ’n’ , 
16 (’q’),’n’, 16 (’q'),’u’) 

120 F0RMATC2X.A1.8X, ’n’ , 16( ’q’ ) , ’n’ , 16( ’q’) , *n* , 

St 16(’q’) , ’n’ ,16(’q’) , *u’) 

DO 130 M=1 , 19 
M1=M 

IFCM.EQ. 5)WRITE(7 , 117) C0(M1:M1), IY34 
IF(M.EQ. 5) GOTO 130 
IFCM.GT. 5) GOTO 122 
GOTO 128 

IFCM.EQ. 10)WRITE(7, 117) C0(M1:M1), IYCEN 
IFCM.EQ. 10) GOTO 130 
IFCM.GT. 10) GOTO 124 
GOTO 128 
Ml*l 

IFCM.EQ . 15) WRITE (7 ,117) C0(M1:M1), IY14 
IFCM.EQ. 15) GOTO 130 
WRITE (7 , 140) C0(M1 : M 1 ) 

CONTINUE 

F0RMAT(2X, A1 ,8X , ’ x’ , 16X, ’x’ ,16X,’x’ , 

St 16X , ’x ’ , 16X , ’x ’ ) 

WRITE (7 , 120) C0(1 : 1) 

WRITE (7, 145) ESC, ESC, ESC 

145 FORMAT C ’ , A1 , ’ *B * , A1 , ’n’ , A1 , ’ [Om’) 

WRITEC7 , 146) XMIN, X14, XCEN, X34, XMAX 

146 FORMAT (8X.F5. 1 , 12X.F5 . 1 ,12X ,F5 . 1 , 12X.F5. 1 , 12X.F5 . 1 
k //37X, ’DIAMETER [microns]’) 

WRITE (7, 161) ESC, ESC 

WRITE (7, 147) ESC, ESC, ESC, ESC, ESC, ESC 

147 FORMATC ’+’ ,A1 , ’ [7;2f 5*. 

St A1 , ’ [8;2f 6’, 

St A1 , ’ [9 ; 2f 7’, 

St Al,’[10;2f 8’. 

St Al, ’ [11 ;2f 9’ , 

St Al , ’ [12 ; 2f :’) 

161 FORMATC’+’ ,A1, ’* C’.Al.’n’) 

DO 168 J = 12,78 
J1 - J/10 
JO ■ J - Jl*10 
IF ( J .LT. 10) Jl“10 
IF( JO .EQ .0) J0=10 
NUMYO « NUMBER (JO: JO) 
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124 


128 

130 

140 



NUMY1 = NUMBER(J1:J1) 


C 

C... HISTOGRAM SCREEN PLOT 
C 

J2 - J-ll 

IF(Y1(J2) .LE.O.O) GOTO 168 
NN2 * 4 

NN1 = IFIX(Y1 ( J2)/YSCAL1) 

II = 2 
10 = 1 

NUMXO = NUMBER(IO : 10) 

NUMX1 * NUMBER(I1 :I1) 

IFCNN1.GT.4) GOTO 162 

IF(NNl.EQ.l) WRITE(7, 163) ESC, (NAME(K) ,K=1 ,7) 

IF(NN1 .EQ.2) WRITE(7, 164) ESC, (NAME(K) ,K=1,7) 

IF(NN1.EQ.3) WRITE(7, 165) ESC, (NAME(K) ,K=1,7) 

GOTO 168 

162 WRITE(7 , 166) ESC, (NAME(K) ,K*1 ,7) 

163 FORMAT( , + * , A1.7A1, * 1*) 

164 FORMATC , A1 ,7A1 , ’2* ) 

165 FORMAT( ’ ♦* ,A1 ,7A1 , ’3* ) 

166 FORMAT (’+’ ,A1,7A1, ’4’) 

DO 167 I • 20, 1,-1 
NN2 - NN2 + 10 
IFCNN2.GT.NN1) GOTO 250 
II = 1/10 

10 = I - 11*10 
IF(I.LT. 10)11-10 
IF (10. EQ. 0)10*10 
NUMXO = NUMBER (10 : 10) 

NUMX1 = NUMBER (II : II) 

WRITE(7 , 169) ESC, (NAME(K) ,K=1 ,7) 

167 CONTINUE 

250 NN2 = NN2 - 10 

NN2 * NN1 - NN2 

11 = 1/10 

10 = I - 11*10 
IF(I.LT. 10)11=10 
IF(IO.EQ. 0)10=10 
NUMXO * NUMBER(IO.-IO) 

NUMX1 = NUMBER (11:11) 

IF(NN2.EQ. 1) WRITE (7, 301) ESC, (NAME(K) ,K=1,7) 

IF(NN2.EQ.2) WRITE(7,302) ESC, (NAME(K) ,K=1,7) 

IF(NN2.EQ.3) WRITE(7,303) ESC, (NAME(K) ,K=1 ,7) 

IF (NN2.EQ.4) WRITE(7,304) ESC, (NAME(K) ,K=1,7) 

IF(NN2.EQ.5) WRITE(7,305) ESC, (NAME(K) ,K=1 ,7) 
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168 

301 

302 

303 

304 

305 

306 

307 

308 

309 
169 
C 

C. . . 
C 


170 
’) 

171 
>F7 . 

172 

= \F7. 

173 

= ’,F7. 

174 

-*.F7. 

175 

= ’ ,F7 . 

176 

= M7) 

177 

178 

179 

198 

199 


IF(NN2,EQ.6) WRITE(7,306) ESC, (NAME(K) ,K-1 ,7) 

IF(NN2.EQ.7) WRITE(7,307) ESC, (NAME(K) ,K=1,7) 

IF(NN2.EQ.8) WRITE(7,308) ESC, (NAME(K) ,K=1 ,7) 

IF(NN2.EQ.9) WRITE(7,309) ESC, (NAME(K) ,K=1 ,7) 

CONTINUE 

FORMAT(’+’ , Al ,7A1 , ’ A* ) 

FORMAT ( ’■*■’ ,A1,7A1,*B') 

FORMAT(’+* ,A1,7A1,»C’) 

FORMAT (’+’ ,A1,7A1, ’D’) 

FORMAT ( ' + ’ , A1 ,7A1 , ’E’ ) 

FORMAT(’+’ ,A1,7A1,’F’) 

FORMATC + ’ .A1.7A1, ’G’) 

FORMAT( ,A1,7A1,»H») 

FORMATC , A1 ,7A1 , *1* ) 

FORMAT( ’ .A1.7A1, * J’) 

OUTPUT MEAN DIAMETERS TO SCREEN 

WRITE (7, 145) ESC, ESC, ESC 
CALL DATE (DAY) 

CALL TIME (TIM) 

WRITE(7 , 170) ESC 

FORMAT ( ’ + ’ , A1 , ’ [3 ;85:f ’ , ’Spatial Distribution 

WRITE(7 ,171) ESC, GCD 

FORMAT( ’ + ’ ,A1 , ’ [5;85:f ’ , ’ Most Probable Dia. 

1) WRITE(7 , 172) ESC, DAV 

FORMAT(’+’ ,A1, ’ [6;85f » ,’ Arithmetic Mean Dia. (DIO) 

1) WRITE(7 , 173) ESC, DS 

FORMAT( ’ + ’ , A1 , ’ [7 ;85f * , * Surface Mean Dia. (D20) 

1) WRITE (7, 174) ESC, DV 

FORMAT(’+’ ,A1, ’ [8;85f * , ’ Volume (Mass) Mean Dia. (D30) 

1) WRITE(7 , 175) ESC, DVS 

FORMAT(’+’ ,A1, * [9;8Sf ’ , ’ Sauter Mean Dia. (D32) 

1) WRITE (7, 176) ESC, NSUM 

FORMAT( ,A1 , ’ [11 ;85f ’ , ’ 

WRITE (7 , 177) ESC, FILE 
FORMAT(’+’ , A1 , ’ [16 ;85f ’ , ’ 

WRITE(7 , 178) ESC, DAY 
FORMAT (’+’ , A1 , ’ [17 ;85f ’ , ' 

WRITE(7 , 179) ESC, TIM 
FORMAT ( ’+’ ,A1 , * [18 ; 85f ’ ,’ 

READ (7, 198) JUNK 
FORMAT(Al) 

WRITE(7 , 199) ESC, ESC, ESC, ESC, ESC 

FORMAT ( IX, Al, ’ [Om’ ,A1, * [H’ ,A1, ’ [?31* .Al, ’ [2J* ,A1 , ’ [?25h’ ) 
CL0SE(2) 


Total Count 


File: 

’ ,A10) 

Date : 

’ ,9A1) 

Time : 

’ ,9A1) 
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CALL SETCMD( ’RUN GRAPH’) 

CALL EXIT 

200 STOP 

END 

C 

C... SUBROUTINE UTILIZED TO SET-UP HISTOGRAM SYMBOLS 

C 

SUBROUTINE CHARGR 

BYTE ESC 

DATA ESC / 27 / 

WRITE(7, 147) ESC, ESC 

147 FORMAT(’+’ ,A1,’P0;33;1;0;0;0{ «???????/GGGGGGG’ , A1 , ’ V ) 

WRITE (7, 148) ESC, ESC 

148 F0RMATO + ’ ,A1 , ’P0;34; 1 ;0;0;0{ «???????/KKKKKKK’ ,A1, ’V) 
WRITE (7, 149) ESC, ESC 

149 FORMAT( ’+’ ,A1 , ’P0;35; 1 ;0;0;0{ «???????/MMMMMMM’ ,A1 , ’\’ ) 

WRITE (7, 150) ESC. ESC 

150 FORMAT ( ’ +’ , A1 , ’P0;36 ; 1 ;0;0;0{ «???????/NNNNNNN’ ,A1 , ’\’ ) 
WRITEC7.151) ESC, ESC 

151 FORMAT(’ + ’ ,A1 , ’P0;37; 1 ;0;0;0{ C /NNNNNNN’ ,A1 , ’\’ ) 

WRITE(7, 152) ESC, ESC 

152 FORMAT ( ’ + ’ , A 1 , ’ PO ; 38 ; 1 ; 0 ; 0 ; 0{ Cooooooo/NNNNNNN’, A1 , ’ \ ’ ) 
WRITE (7, 153) ESC, ESC 

153 FORMAT ( ’ + ’ , A 1 , ’ PO ; 39 ; 1 ; 0 ; 0 ; 0{ Cwwwwwww/NNNNNNN ’ ,A1 , ’ V ) 
WRITE (7, 154) ESC, ESC 

154 FORMAT ( ’ + ’ , A 1 , ’ PO ; 40 ; 1 ; 0 ; 0 ; 0{ «{{{{{{{/NNNNNNN ’ ,A1, ’\’) 
WRITE (7, 155) ESC, ESC 

155 FORMAT( ’ + ’ ,A1 , ’P0;41 ; 1 ;0;0;0{ C>»»»/NNNNNNN’ ,A1 , ’ \’ ) 
WRITE(7 , 156) ESC, ESC 

156 FORMAT(’ + ’ ,A1,’P0;42;1;0;0;0{ C /NNNNNNN’ ,A1 , ’V ) 

WRITE (7, 157) ESC, ESC 

157 FORMAT( ’+’ , A1 , ’PO; 17 ; 1 ;0;0 ;0{ CWWWWWWW/ ??????? ’ , A1 , ’ V ) 
WRITE (7 , 158) ESC, ESC 

158 FORMAT (’+’ , A 1 , ’ PO ; 18 ; 1 ; 0 ; 0 ; 0{ «[[[[[[[/???????’ ,A1 ,’ \’ ) 
WRITE (7, 159) ESC, ESC 

159 FORMAT ( ’+’ ,A1,’PO;19;1;0;0;0{ «]]]]]]]/???????’ ,A1 , ’\’ ) 
WRITE(7 , 160) ESC, ESC 

160 FORMAT ( ’ + ’ , A 1 , ’ PO ; 20 ; 1 ; 0 ; 0 ; 0{ « /???????’ ,A1 ,’ \’ ) 

WRITE(7 , 161) ESC, ESC 

161 FORMAT (’+' ,A1,’P0;21;1;0;0;0{ CNOOM««]/??????? ’ , A1 , ’ \ ’ ) 
WRITE (7, 162) ESC, ESC 

162 FORMAT ( ’ + ’ , A 1 , ’ PO ; 22 ; 1 ; 0 ; 0 ; 0{ «*CCCCCC/??????? ’ , A1 , ’ V ) 
WRITE (7, 163) ESC, ESC 

163 FORMAT( ’+’ , A1 , ’P0;23 ; 1 ; 0;0;0{ CPXTTTRP/ ??????? ’ , A1 , ’ \ ’ ) 
WRITE(7 , 164) ESC, ESC 

164 FORMAT ( ’ + ’ , A 1 , ’ PO ; 24 ; 1 ; 0 ; 0 ; 0{ «PPPPPPM/ ??????? ’ , A1 , ’ \’ ) 
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165 


166 


WRITE(7 , 165) ESC, ESC 

FORMAT( , +’ ,A1,’P0;25;1;0;0;0{ CMPPPPPM/???????’ ,A1 . 'V ) 
WRITE(7 , 166) ESC, ESC 

FORMAT( , A1 , ’PO ; 26 ; 1 ;0; 0 ; 0{ CMPOOOPM/??????? ’ ,A1 , ’ \ ’ ) 
CALL EXIT 
STOP 
END 
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Section 10 


APPENDIX C.3: MOD-1 Nozzle Input Pressure Determination 


C 

SUBROUTINE PRESSURE(NCHCK) 

C+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

+++++++++++ c 

C SUBROUTINE TO DETERMINE WATER AND AIR PRESSURE FROM C 

THE OMEGA TRANSDUCERS USING THE AXV11-C A-D BOARD. C 
C+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
♦+♦+♦+♦++♦+ L0GICAL*1 I KEY 

BYTE ESC 
DATA ESC / 27 / 

CALL IP0KE( "44 , "10000 . OR . IPEEK ("44) ) 

C 

C. . . THE FOLLOWING ARE THE APPROPRIATE OCTAL VALUES TO 
C. . . BE STORED IN THE CSRs OF CH. 1 A 2 TO START AN 
C. . . A TO D CONVERSION OF THE PRESSURE TRANSDUCERS. 

C 

ISTRT1 = "415 
ISTRT2 = "1015 
C 

IF(NCHCK.Eq.l) GOTO 12 
WRITE(7, 10) ESC, ESC, ESC 

10 FORMAT( ’ ,A1,’ [2J* , A1 , ’ [?251 ’ , A1 , ’ [H ’ ) 

WRITE(7, 11) ESC 

11 F0RMAT(1X,A1 , ’#6’ , ’ PRESS "C" to continue’) 

C 

C... 17770400 IS THE CSR (CONTROL STATUS REGISTER) FOR CH. 1 

C 

12 CALL IPOKE(" 17770400 , ISTRT1) 

13 ICHK » IPEEK ("17770400) 

IF(ICHK.NE."614) GOTO 13 

C 

C... 17770402 IS THE DBR (DATA BUFFER REGISTER) FOR CH. 1 

C... AND IPW IS THE WATER PRESSURE (DIGITAL VOLTS). 

C 

IPW = IPEEK ("177770402) 
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14 


C 

C. . - 
C 


18 

t 13,’ 

19 

ft 13, 

20 

21 

25 


CALL IPOKEC" 17770400, ISTRT2) 

ICHK * IPEEKC" 17770400) 

IF(ICHK.NE. "1214) GOTO 14 
IPA * IPEEKC "17770402) 

IPA = -1.189 + 0 .4S98*IPA 
IPW = -1.95 + 0.4598+IPW 

OUTPUT PRESSURE VALUES TO TERMINAL SCREEN 


WRITE (7, 18) ESC, ESC, IPW. IPA 
WRITEC7.19) ESC, ESC, IPW, IPA 
IF(NCHCK.Eq.l) GOTO 25 
IKEY = ITTINRO 
IFCIKEY.EQ. ’C’) GOTO 20 
GOTO 12 

FORMAT ( ’ + ’ , A 1 , ’ [22 ; If ’ , A1 , ’ #3 * , ’ WATER PRESSURE = 
AIR PRESSURE = ’ ,13) 

FORMAT ( ’+* ,A1, ’ [23; If ’ ,A1, *#4* , * #WATER PRESSURE 
AIR PRESSURE = ’,13) 

CALL IPOKE ( "44 , "167777 . AND . IPEEK ( "44) ) 

WRITEC7 ,21) ESC, ESC 

FORMAT C * + ’ ,A1 , 1 CH’ ,A1 , ’ [2J’ ) 


RETURN 

END 
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Section 11 


APPENDIX C.4: PSP Magnification Correction Factor Determination 


+++++++♦ c 

C FUNCTIONS TO DETERMINE CORRECTION FACTORS FOR 

C MICRON TO PIXEL FACTORS WHICH DEPEND ON X AND Y. 

C 

+♦♦+++++ FUNCTION XFACT(MAG,XPOS ,YP0S) 

IF(MAG.EQ.SOO) GOTO 100 
C 

C... HIGH MAGNIFICATION X-CORRECTION 

C 

IF (XPOS . LE . 100 . 0) XFACT=0 . 977+YP0S+8 . 09E-05 


IF (XPOS . GT . 100 . 0 . AND . XPOS . LE . 150 . 0) 

XFACT=0 . 974+YP0S*2 . 60E-05 

IF(XPOS . GT. 150 . 0 . AND . XPOS . LE . 200 . 0) XFACT=0 . 967-YP0S*8 . 12E-07 

I F (XPOS . GT . 200 . 0 . AND . XPOS . LE . 250 . 0) XFACT=0 . 96 1+YP0S*4 . 73E-06 
IF (XPOS . GT . 250 . 0 . AND . XPOS . LE . 300 . 0) XFACT=0 . 961- 
YPOS+5 . 46E-05 IF (XPOS . GT . 300 . 0 . AND . XPOS . LE . 350 . 0) 

XFACT=0 . 948-YP0S*3 . 72E-05 

IF (XPOS . GT. 350 . 0 . AND . XPOS . LE . 400 . 0) XFACT=0 . 943-YP0S*6 . 80E-05 
IF (XPOS . GT . 400 . 0) XFACT=0 . 920- YP0S*2 . 58E-05 
RETURN 

100 IF(XP0S.LE. 100.0) XFACT=2 . 21+YP0S*0 . 290E-04 

C 

C... LOW MAGNIFICATION X-CORRECTION 

C 

IF (XPOS . GT . 100 . 0 . AND . XPOS . LE . 150 . 0) 

XFACT=2 . 22+YP0S*0 . 803E-04 

IF(XPOS . GT. 150 . 0 . AND . XPOS . LE . 200 .0) XFACT=2 . 16+YP0S*O . 679E-04 

IF (XPOS . GT . 200 . 0 . AND . XPOS . LE. 250 . 0) XFACT=2 . 16+YP0S*0 . 442E-07 
IF (XPOS . GT . 250 . 0 . AND . XPOS . LE . 300 . 0) XFACT=2 . 16- 
YP0S*0 . 947E-04 IF (XPOS . GT . 300 . 0 . AND . XPOS . LE . 350 . 0) 

XFACT=2. 11-YPOS*0.306E-07 

IF (XPOS. GT. 350.0. AND. XPOS. LE. 400.0) XFACT=2 . 10-YPOS*0 . 124E-03 
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IF CXPOS. GT. 400.0) XFACT=2 .07-YPOS*0. 135E 03 
RETURN 
END 

FUNCTION YFACTCMAG,XPOS,YPDS) 

IF (MAG. EQ. 500) GOTO 100 


C 

C. . . 
C 


HIGH MAGNIFICATION Y-CORRECTION 


IFCXPOS . LE . 100 . 0) YFACT=0 . 977-YP0S*9 . 17E-05 
IF (XPOS .GT. 100 . 0 . AND . XPOS .LE. 150.0) YFACT-0.981 
YPOS* 1 . 24E-04 IFCXPOS . GT . 150 . 0 . AND . XPOS . LE . 200 . 0) 

IFCXPOS GT^OmNrXPOS.LE .250.0) YFACT=0.990-YP0S*1.63E-04 

i?apSs G? 2sJ S.AKD.XP0S.LE.300.0) YF.CT-1.000-YPOS.1.96E-04 

IF CXPOS. GT. 300.0. AND. XPOS. LE. 350.0) YFACT=1.014 

YP0S*2 . 19E-04 IF(XP0S - GT - 35 °-°- AND ‘ XP Sn L ni 4 Y?A?T-l 029- 

YFACT=1 . 027-YPOS*2 . 63E-04 IF CXPOS. GT. 400.0) YFA 

YP0S*2 . 69E-04 
RF.TURN 


100 

C 

C. . 
C 


IFCXPOS . LE . 100 . 0) YFACT-2 . 10-YPOS*0 . 183E-03 
LOW MAGNIFICATION Y-CORRECTION 


IF CXPOS. GT. 100.0. AND. XPOS. LE. 150.0) YFACT-2. 15- 

YP0S*0 . 268E-03 IFCXPOS . GT . 150 . 0 . AND . XPOS . LE . 200 . 0) 

YFACT=2 12-YP0S*O . 315E-03 

IFCXPOs'gT. 200.0. AND. XPOS. LE.2S0.0) YFACT=2. 13-YP0S^.313E-03 
( IF CXPOS. GT. 250.0. AND. XPOS. LE. 300.0) YFACT=2 . 15-YP0S*0 . 397E 
IF CXPOS. GT. 300.0. AND. XPOS. LE. 350.0) YFACT-2. 18- 
YPOS.0.464E-03 ^ (XPOS . CT . 350 . ° ■ AND . XPOS . LE .400^) 

YFACT=2,19-YP0S.O,6O5E-O3 IF(XPOS.GT.400.0) YFACT-2. 18 


YP0S*0 . 509E-03 
RETURN 
END 
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Section 12 


APPENDIX D: Mean Diameter Calculations 


Arithmetic Mean Diameter (AMD) 


0(10) = 




E £ 


1 


Area Mean Diameter (ArMD) 


D( 20) = 




EL ^5 

Eili «. 


Volume Mean Diameter (VMD) 


(12.1) 


£>(30) = 


N 


E El M? 

2-/t=i n * 


Sauter Mean Diameter (SMD) 


E£i M? 
E£i M? 

where 

N = total number of bins 

ni = counts per bin 

di - diameter for size class i 



( 12 . 2 ) 


(12.3) 
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Section 13 


APPENDIX E: Cole- Palmer Flowmeter Calibration Data 


Scale Reading 
30 
40 
50 
60 
70 
80 
90 
100 
110 
120 
130 
140 


Flow-rate (gpm) 1 
0.022 
0.031 
0.040 
0.050 
0.059 
0.069 
0.078 
0.090 
0.100 
0.110 
0.120 
0.130 


♦Calibration values were verified by replication. 
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Section 14 


APPENDIX F: OMEGA Pressure Transducer Calibration Data 


Table 14.1: S/N: 850502 
Standard Pressure Output Voltage 


Corrected to psia 

(millivolts) 

0 

9.20 

10 

15.00 

20 

21.70 

30 

28.30 

40 

34.80 

50 

41.45 

60 

47.95 

70 

55.20 

80 

61.70 

90 

68.60 

100 

75.20 

110 

80.70 


Table 14.2: OMEGA S/N: 850311 
Standard Pressure Output Voltage 


Corrected to psia 

(millivolts) 

0 

9.80 

10 

16.00 

20 

22.90 

30 

29.10 

40 

35.75 

50 

42.60 

60 

49.10 

70 

56.20 

80 

62.60 

90 

69.40 

100 

75.80 

no 

81.50 


* Calibration values were verified by replication. 
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Output Voltage (mV) 


S/N 850502 



Standard Pressure (pela) 


Figure 14.1: OMEGA Pressure Transducer Calibration 
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Output Voltage (mV) 


S/N 850311 



Figure 14.2: OMEGA Pressure Transducer Calibration 
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